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Vol. XII, 1961 


On the Theory of Hot-Water Drives 


By Cuarrrs L. Carney, Jr.1) and Ropert P. GitBerr?), 
Pittsburgh and East Lansing, USA 


1. Introduction 


In this paper we are going to investigate a hot water drive in an oil-bearing 
tock formation. The hot-water drive is one of several thermal recovery methods 
[1], [2], [8]*) used by the petroleum industry to produce the residual oil in a 
reservoir. (The essential phenomenon all such methods are based on is that by 
heating the reservoir the viscosity of the oil is lowered and may be more easily 
produced.) A disadvantage of this process is that initially cold oil must be 
displaced, unless of course there is ample stratification or pore-space in the 
formation for the hot water to flow through. 

There have been a number of theoretical hot-water drive studies made 
previous to this one. Most assume a piston like displacement by the encroaching 
water; however, they differ on the modes of heat transfer. SCHUMANN [3], 
KLINKENBERG [4], and WALTER [5] have studied a model in which the heat is 
transfered by 

(1) convection of hot liquids, and 
(Model A) (2) heat exchange between the solid matrix and liquids which are 
at different temperatures. 


JENKINS and ARONoFSKY [6] obtained a solution to the case where heat is 
transfered by 
(Model B) (1) convection of hot liquids, and 
(2) conduction in the solid matrix. 


They also compared both models against éxperimental data, and concluded 
that the mechanism involved was in a large part similar to the heat transfer of 
model B; however, in the case of inhomogeneous sands ‘field experience would 
be needed to determine whether model B or a combination of both models 
would sufficiently describe the transfer mechanism’. 

In this paper a solution to the combined models (4 plus B) is sought. This 
is done by performing an iteration on the solution to model 4; this solution 


1) Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania. 
2) Department of Mathematics, Michigan State University, East Lansing, Michigan. 
3) Numbers in brackets refer to references, page 496. 
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will be more close to the actual situation than model A and in the case o 
stratified formations should be more accurate than the solution to model B. 
Nevertheless, greater accuracy always may be obtained by continuing the_ 
iteration process. 


2. Mathematical Formulation of the Problem 


It would be well to list here the symbols used.in the following work. 


Variables 


x distance measured from the input well, 
t time measured from the instance when flooding begins, ; 
T, temperature in the solid, ; 
T,, temperature of the water. ; 


Constants 


v_ velocity at which a piston-type water front moves, 
h,, heat transfer coefficient from water to solid, 
h, heat transfer coefficient from solid to water, 
K, conductivity in the solid, 
k, diffusivity in the solid = K,/(1 — f) Pc, , 
P. density of solid, 
c, heat capacity of solid, 
f porosity, 
T, temperature of input water. 
y=h, = Se ee (ie ak b= orrey, a= nie 


The partial differential equations for heat transfer by conduction in the 
solid, convection in the liquid, and a skin-effect between the solid and liquid are: 


OT, OT, 
rae sir) re sy —h, (Ty <<, T,) 7 | 
and : 
ay, rT, 
Seb, Se mh, (Py), | 
where 
Me YG 


J (1 Ss f) ey es 
is the thermal diffusivity of the solid. 
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oT, pT, 

ot dx bw (Ty — Th) , | 

ged ‘ie (2) 
SS ae h, ee T,) = | 


‘T,, and T, both satisfy the hyperbolic partial differential equation 
§ 
(Cl o?T oT Cia 


pt ee ee Vr He) a 0 (3) 


which has the characteristics x = constant, ¢ — x/v = constant. Equation (3) 
tay be reduced to normal form by the transformations 


: Semele e 
and (4) 


Equation (2) then becomes 


Ou Ou 
i = u, ‘ A = Uy A (5) 


which in turn satisfy the standard equation 


0*u 


oy rab giana (6) 


We may solve this equation readily by means of the Laplace transform method. 
Transforming with respect to the variable z, one obtains 
Ou = Ou oe 
by s) — uly, s) + ree 0)=0. (7) 


If one uses the boundary conditions 


oT 
T(y, 0) = 5°-(y, 0) =0 
at the leading-edge (z = 0), and 
T(0, z) = T, 


at the input well (y= 0) then Equation (7) has the solution 


ii(y, s) = 22 ev" (8) 


Ss — 


“ 
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ite iat antes 


The inverse transform of u(y, s) is 
2Vyz ‘ 


uly, 2) = Te + To f Iyiy) 7" dy, (9) 


0 | 
where J,(y) is a modified Bessel function of the first kind and first order. | 
Equation (9) may be obtained from the tables [7] by use of a convolution. 

Returning to Equation (1) and using the usual transformation to normal — 
form one obtains : 


OUy 
“spe ae 

and (1) 
Ous Re doe ( Ou, Ou, ) 
0z ae hve \ Ope oy boat k 


It may be seen that w,, and u, satisfy the partial differential equation 


07u 
Oy 02 


(10) 


ane es (55 0?u | 


oy? oy? oy 


where a! = k, h2,/h,v®. This equation may be reduced to an ordinary differen- 
tial equation in y by taking the Laplace transform and using the above men- 
tioned boundary conditions. This subsidiary equation then may be written as 


Ou 


0 8) 254.5) + (1 —sa) 4,5) -aly,s)=0. (A) 


By introducing the change of variables 
- Ou 0? 
wi() =H, 8), 0) = 5°0.5), w= F508). 12) 


Equation (11) is reduced to the following system of ordinary, linear-differential 
equations 


wi = f(y; Wy, We, Ws) = We, | 
Wy = falV; Wy, Wy, Ws) = W3, (13) 
W, = fa(y; Wy, Wa, Ws) = 2 ws + (sa — 1) we— wy. | 


A unique, continuous solution exists for this system equations which passes 
through the point (y®, wf, wg, wg), if there is a region containing this point in 
which the /x(y; w,, We, ws) (K = 1, 2, 3) are continuous, and the following 
requirements hold: 


[fa (ys OF, we, WE) — f(y; Wy, We, Ws) | 


| (13a) 


SL (|wt — w, | + | wf — w,| + | wf — w, |) 


f 
“id 
“4 
f . 
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where L is the Lipschitz constant. 


} 


“ | fe(Vs @1, Wy, Ws) | SM for (y, w, w,,w;)ER (13b) 


: (that is the fx are bounded in this region). 
_ The region R may be bounded by the inequalities 


ly=y| Sa, and fw ve | <5, <0. (13c) 


We notice that these three conditions are obviously satisfied from Equation (13), 
and proceed to define the functions 


¥. 
W044) = el) + feos er ao» Wes) AY (14) 
ye 
where K = 1, 2, 3, and v= 1, 2,3,.... It follows from the Cauchy-Lipschitz 
theorem that these successive approximations converge to the solution of our 
system of equation, and in particular that the function 


w4(y) = uly, s) = lim @,, ,(9) (15) 


satisfies Equation (11)*). We proceed by computing the first iteration to w,(y) 

‘starting with w{(y) = T, e¥!’/(1 — s), the solution to Equation (7). Obtaining 
the first iteration to w,(y) from the system of Equation (13) is equivalent to 
solving 


y/s 
(a) wy —2u,+ (1+sa) W,=—-——, (16) 
and 
y 
W, = [w, dy . (16a) 
Vo 

The solution to the homogeneous equation associated with (16a) is 

wy = A elt tVsa)y 4 B elt-Vsa)y. (17) 


4) For s fixed and real one may visualize a finite w — y region as contained in a compact 
subset of a w—y ‘plane’, hence by the Heine-Borel property it may be covered with a finite 
number of rectangles. By starting in an initial rectangle and applying the Cauchy-Lipschitz criteria 
one may then proceed from rectangle to rectangle and obtain a continuous solution w,(y) = u(y, s) 
for all finite values of y. However, this does not ensure that our limit function w,(y) = a) 


is a Laplace-Transform. For instance, let us consider the sequence of functions {@1, (y) = u(y, 8) hy 

and let us suppose that for finite v each %,(y, s) is analytic for Re s > c,, and that consequently an 

inverse Laplace Transform may be obtained by a Bromwich contour integration. Here, each 

function %,(y, s) which .corresponds to a finite y, may be transformed into a function U,(y, 2). 

There is the possibility, however, that the Eats —> oo; in this case it is not self-evident that 
> 


u(y, s) =lim %»(y, s) is a Laplace transform, or that it corresponds to the function U(y, z) = lim U,(y,z). 
v—>0O v—>0o 


This question is not discussed in this paper, but will be investigated in a further note, 
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by variation of parameters the solution to the linear equation may be found 


that is is 
Wee R= Lee 
we(v) = cy gli+Vsa)y 4 Co gi-Vsa)y _ Cie (18) 


The constants c,, c, are determined from the boundary conditions, at the leading- | 
edge (w, = 0, y =/,, t) and the input-well . 
i | 
(w= m=, y=), | 
| 
to be [ 
L, [2sae7Vsahwt 4 pt/s (4 — ¢-Vsahwt)| t 
y= —_ ’ 
* 45a (1—s) sinh(//s « h,, ) | 
and ] 
tie —T, [2s aeVsahwt 4 pils ( 1 — ¢-Vsahwt)] : (19) 
: 45a (1—s) sinh(/s « h,, ) 
Ou, ee oT ner: 
an ey (T+ a ee | 


at v0, for 2 >-b. 
By rearranging terms in (18) and (19) one has 


ou sinh [(y — hy #) Vor s?”] 
(Vy, s)=w = =f) eY : ees SS 
oy (y ) 2(¥) b le (s — 1) sinh iy , Vo sil2] | 
4 ella cosh[(y — Ayre 2) Ya s/2] _T, etaty (20) 
2sa(s — 1) pene st a aw 
2 


The inverse Laplace transform may be found by use of the tables [7] and the 
convolution theorem; one obtains 


Ou kh, he 
one = Yore | sw 
By 02 2) = Tye? EB 


“EVE [lente VEl: 
ot V h, of “4\20e Va, | howe) 
0 


nH 


4 be « : ”.)| dn at 


es 


0 £ [he 
a ea m1 (f LE 


gn€ 


= eae 
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; ky hy ie ad : Rk, hy 2 
where 
. 2 
0,(v | T) = (— iz) S(—)n e7mv-W2+n) |e ‘ 


n= —O0O 


(22) 


6,(v'|-z) = O,(8 | t) =(- 1) ABS" i O-12 +7 fp 


n=—0o 
are theta functions. We recall that y, z are the independent variables and ¢ 
appears just as a parameter because of the boundary conditions, hence one may 
write 


2 z 
u(y, z) = f(z) + T, ert? oa f e~* Ig(2 Cu?) dl 
Z 0 


Zz 


fl kg F YS Ny t in Chk; =e, 
ee on en)” 
0 
y —hygt/2 | 4 aik,(€ — 7) ee 
bg fe nan | SE ana 


oe (n—£) —hyt inch, 
+ Ti we fd rer acc r) ao dy 


(€—¢) 1/2 
ai {fe T9(2 7”) 


E—hyt/2 | 4 ath, (C — 0) 
x 8, ( Tay € 2 h,(v #2 ) dn ag ag 


where f(z) must be determined by the y boundary conditions. The two solutions 
u,(y, 2) and u,,(y, 2) may be distinguished at this point by using the conditions 


u,(0, 2) = T,,(0, 2) € = T, e (Yo= 9), 
and 


tt(h, t, 2) = T, ryt, 2) = 0 (Yo = hut) - (24) 


By choosing y, differently is each case, /,(z) and f,,(z) may be calculated from 
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the equations 


fyol2) = & T, — & Ty 


‘ 0, ( 1 i re, DY dn at 


fs) = —1, eh" | 24, [et12o% 4 
2 AG 


te J a(0 


Zz 


Wee 25b 
Tan oe eee ke | | e* Io(2 9") a 
0 0 


x 6, (| sa) an a]. 


Equations (21) and (25) represent a first approximation to a solution for 
the hot-water drive problem. 
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On the Instability of Taylor and Dean Flows under the 
Influence of a Circular Magnetic Field 


By Rosert H. DYER and Davip E. Epmunps, Cardiff, Wales?) 


1. Introduction 


In a recent paper, YIH [1]?) has shown that ring vortices may be produced 
in a stationary fluid confined between two concentric circular cylinders when 
an electric current of suitable constant density passes axially through the fluid, 
and a line current passes along the axis of the cylinders. The instability is 
brought about by the centripetal electromagnetic body force. Cases for which 
the production of similar vortices is of a dynamical character are, of course, 
well known (cf. TAYLOR [2], DEAN [3)). In [1] sufficient conditions for stability 
are obtained by integral methods similar to those of SYNGE [4], and detailed 
results are given when the spacing between the cylinders is small compared 
with the radii. 

In this note an extension of Yru’s work is made to consider the problems 
arising when the undisturbed fluid is not stationary but has a basic velocity 
around the cylinders. This may be achieved either by rotation of the cylinders 
or by the application of a pressure gradient round the cylinders. In such 
problems two sources of instability are present, and a discussion of their joint 
effect is given here. Y1H in his paper notes that a simple experiment to demon- 
strate the phenomenon does not appear possible owing to the great magnitude 
of the current that would be required. The position is slightly better with 
regard to the present work in that we should require only a current large 
enough to produce instability at a TAYLOR number perceptibly different from 
the usual non-magnetic one. Unfortunately this current would still be formida- 
bly big, if mercury were used as the fluid. 


2. The equations of the problem 


We shall consider the motion of a viscous, incompressible, electrically 
conducting fluid in the space between two infinitely long coaxial circular 
cylinders, whose equations are r= Ry, 7 = R, (Rz > R,) referred to a system 
of cylindrical polar coordinates (7, 0, z). The equations governing the motion 


1) Department of Mathematics, University College, Cardiff, Wales. 
2) Numbers in brackets refer to the references, p. 507. 
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of such a fluid are, in component form 


2B 8) -— erm 2B + EB) 
a( Bas ttt) 15 x0r(rme t+ 3 95) 
$e pay 


DHg» up H, Dug cm si0n +n (V2 H, Ho a. 2 ey, 


Dt + Y Dt y2 v2 06 


DH, Du, 2H 
Des ti a ae 
Ou, Uy 1 0ug , Ou, 7 
Or mts Y a 06 + 0z Oe ”) 
OH, iE 1 0H¢ Oe 8 
Paani Ma aS aT Pee ®) 
where 
= w|i? 
= p sir 82 + Q 2 - 
KRY) 0 | uw 2 0 
a Ot oie Or te yy 00 te On, 
Duk. 0 Figs Die 0 
Df ne 0 ago a0 coer ‘ 
0? Loo Layer 0? 
ee se. ee ae 
lee or? Y Ov a v2 062 Oz? ’ 
and 


N= (490609) 5+. 


In these equations, ¢ is the time, @ the density, « the permeability, v the kine- 
matic viscosity, o) the conductivity, and 2 the gravitational potential per unit 
mass; U,, Up, ¥, are the components of the velocity u and H,, Hg, H, are the 
components of the magnetic field H in the radial, transverse and z-directions 
respectively. The current density j is such that 


j=o,(E+pvuxaA), corlH=42j. (9) 
E being the electric field. 
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It is easy to verify that equations (1) to (9) have the following two station- 
ary solutions: 


(a) b= 64,209 Ug=V =agr+b,7', (10) 
: We vaked ee 
Pant = Oe Hy=H=22(+— + fo" (11) 


where 
Q, — Q,) Ri RR 
Ri — R? 


Eoagt| 
. a 
[= J =n Bin. 


The velocity field is produced by rotating the inner and outer cylinders 
with angular velocities 2,, Q, respectively; the magnetic field is the result of 
passing a current J along the centre line of the cylinders, i.e. the z-axis (J is 
taken as positive if it is in the positive z-direction), and a constant current of 
density j, between the cylinders in the z-direction. There is no interaction 
between the velocity and magnetic fields. 


(b) y= z= 0, My =V=K(C(e me log | 12) 


where 
dee -1 R > Rage? 
(aay log R,’ a ne 
The magnetic field is taken to be as in (a), while the velocity field is generated 
by a constant pressure gradient in the 6-direction. 

The stability of these two types of flow will now be investigated, and the 
two problems thus arising will be referred to in what follows as the TAYLOR 
problem and the DEAN problem respectively, because of the analogies with the 
work of TAYLOR and DEAN in the absence of a magnetic field. We shall establish 

‘the basic stability equations without actually specifying V, so that up to this 
point the work will apply to both problems. 

Let the stationary flow given by (11) and either (10) or (12) be slightly 
perturbed, so that the components of the velocity and of the magnetic field of 


the disturbed motion are 
Oop Via tes Ces h,,H +hg,h,, (13) 
respectively. The perturbation in the magnetohydrodynamic pressure ¥ will be 


denoted by yz’. Assuming that the disturbance is independent of 0, we find on 
substituting (13) into equations (1) to (6) and linearising: 


H 1p.0x. 1 
OV, = 2. VA V6 _ ERES he ry 0. = + y (V2 _ =z) Us ) (14) 
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Ove _ IV ill yt dH add ees sie 
ees, sera Bayi ez cier  calld ps) % ( : 
Oy Aor 2 16 
ot 0z a PVE Mes ( . 
Np _ ret 17\4 
Ot n(V ya) Ie . 
bheceut dH Ses 6 | 
Sa ee et ae ee ee (Vr p) te a | 
Oh, 2 1 q 
where now 
¥ 0 I RRya, 02 
ove a 
= or? 7 Y OY a Oz? 


The form of equations (17) and (19) leads us to take both A, and h, to be 
zero functions. Equation (7) enables a stream function y, to be introduced, 


so that ere 
0p i ee 


¥ 
and we specify the type of perturbation to be considered by writing 
Yi = —tp(r) exp(ot+iAz), vg=d(r) exp(ot+7Az), 
ha = h(r) exp(ot+7iA2). 


Elimination of y’ between (14) and (16) then gives 


(L—2— 2) (L—-ao=- 22 [va-#E Al, (20) 

while (15) and (18) become 
(LB 2)6-—1p (2 4 7), (21) 
n(L — 2 — <) here (= xe =) ; (22) 


where 


The velocity boundary conditions are provided by the usual no slip require- 
ment, and are 


Nother dy 
6=9=% <0 at r=, 1HR. (23) 


If we suppose that the material of which the cylinders are made isan insulator, 
then 7,=0 on the walls. Taken together with (9) this gives the magnetic 
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boundary condition as 


Our problem is thus to solve equations (20), (21), (22) subject to the boundary 
conditions (23) and (24). 


3. The Taylor problem 
_ Here the basic velocity field is specified by (10). Defining dimensionless 
quantities 
Ye ee Ta Aig 0 ee yee ay 


i Ps 
6 R, Rig) ie, ta 
Te gg or tak (ae ame 


and then dropping the primes on 7’ and o’, we find that (20), (21), and (22) 
become 


b B 
(Lae = al lL #)y=—2k(a+—)ot(1+a)h, (25) 
(L—k?-—o)v=—2aky, (26) 
(E-#- 2”) nan », (27) 
where 
R (B a? — 1) b R a? (B — 1) 
rs a 1” a2?—1 ’ 
Q, Ri _ _ As 
R= 11, pag «= Re (28) 
— #2 Rt code ite 248 
ced hel’, n 


Proceeding exactly as in [1] and [4] it is possible to obtain from (25) to (27) 
the relation 


REN (Io + 2k 1, + hI + Jot Js) + Hot A, 
+0,{8N BIp+ t+ 1) +> Hf =0 


in which 

ra g a 

Iy= f rly Par, = | er cae a | dr, I,= "| r|Ly|dr, 
1 : 1 ul 


, b 
“dr, T= felts) lyre. 


F 
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a 
ee B qe 4 i 3 2 
Hy = {r(i+ 3) | Lig iy dre, Hs fr (1+ 4) [a Par 
£ 
and a, is the real part of o. All we can deduce from this, however, is that the 


motion is stable if both 
(i) OPE cet 8 ON Pee ON 


and 
(Gof zine], Ry er ef Ss eg (ey) ' 


hold, i.e. if both the Yih and Synge conditions hold. | 
We now turn to the solution of equations (25) to (27). These are too difficult 
to solve as they stand, but when the spacing d = R, — R, between the cylinders - 
is small compared with the radii the problem simplifies. If the cylinders are 
rotating in the same direction we may approximate a + b/r? = w in (25) by 
some suitable value w such as 
7 ? = ae 
Oa aes fo cs (1+ aa 
and similarly if the two applied currents are in the same direction, i.e.if -B<1, 
then the quantities 1 + B/r? and 1/r? occurring in (25) and (27) respectively 
may be approximated by 1+ B/r§ and 1/r%. Defining new dimensionless- 
parameters 
ee ee tee RY ak 
o= 7 (y=), m=k R,’ a= 0 (2) 
and neglecting terms of order d/R,, equations (25) to (27) become, under the 
above approximations, and with D = d/dé, 


(D? — m? — o,) (D? — m*)\ p= (=) & 2m (=) wov+ (1 + = =) n|, 


Ry 


ad 
D* — m? — o,) v= —2am(—)y, 
( 1) (a) Y 
(D? — m? — 0, »/n) h = “> ke 


Assuming the principle of exchange of stabilities, we thus find that 
(D2 — m2) y = — m2? (T+ Q) yp 
where T = — 4w a (d/R,)‘ is the Taylor parameter, and 
ie. (1 + B/ro)N(d/R,)4 


2 
‘9 


The boundary conditions are 


p= Dy= (DP —wAp=0 at f= 077, 
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and so the problem is equivalent to that of PELLEw and SouTHWELL [5]. Using 
their solution it follows that instability will first set in when 


f+ O=1707°8. 


Hence the sign of Q determines whether the magnetic field stabilises or desta- 
bilises the flow. 

In [1] it is shown that to produce vortices in mercury at rest between two 
cylinders with 

B=-<, d=1cm, Rk,=10cm 

the total current through the annular space would have to be 1-62 x 10° amp. 
Our result shows that to depress the critical TAYLOR number T to a value 50 
below its non-magnetic value, the total annular current required would be 
roughly one sixth of the current mentioned above. Thus for the predicted effect 
to be detectable a very large current would still be required, if mercury were 
the fluid used. As pointed out by Yiu, this current would be greatly reduced 
if it were possible to replace mercury by a highly ionised gas with small density 
and high conductivity. 


4. The Dean problem 


Here we have to solve equations (20) to (22) with the basic velocity field 
defined by (12). As in § 3 we shall consider only the case where 


J) 
Ce Bey ee (R, + Rs) 


and then we have approximately 


, 3 
V=5 V, (1-424 
with : 
Se ity 1 Piece, op: 
ea Tamer Aa 120vR, 00° 
Writing 
A m - y ets 2a Br. 
t= V_v, 6 er ee Gant ae === ee -foah 


and neglecting terms of order d/R,, we find that the perturbation equations 
become (assuming the principle of exchange of stabilities) 


(D2 — m2)? y= mT (1-4 x?) v —m* Th, (29) 
(D2 —m?)v=—2xy, (30) 
(D?—m)h=y, (31) 
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where B and N are defined in (28), and 
18 Rd Vy 4 a d\4 tee. 
7 R= 25, --Ni+B\(p). D=—=. 


The boundary conditions are 


y= Dy=(D?—m)yp=0 at x=+41/2. 


Equations (29)-(31) are now solved using the method of REIp [6]. _ 


assume that y can be expressed in the form 
y PD she C,(x) + B, S,(x)] (32) 
n=1 ; 
in which ; 
C _ coshé,#% COSA, % (x) = sinh fy, * SIN My ¥ 
JC epranry ay Tae Mk an sinh p,/2 sin ,/2 


where 4, and yw, are the positive roots of the equations 


ih a i) 1 
tanh 7 A+ tan>4=0, coth > m—cot > m=0. 


The functions C,,and S, of course satisfy the boundary conditions on y. Substi-. 
tuting (32) in (30) and (31) we find solutions of the resulting equations for v 
andh which satisfy the boundary conditions. If we substitute these expressions 
for v and fA, together with (32), into (29),and multiply the resulting expression 
in turn by C,(x) and S,(x) and then integrate between x = — 1/2 and x =1/2, 
we obtain the following system of equations for the A, and the B,: 


y q=1,2,...) (33) 
1 
2 (A, ng) pe ee B,,) 0 | 

where 
Ys 

Ang = (By + 1) Bug — 20 Xyq t mT, [fy Palm) ++ a saa (OH Dag + Xnq)| 
1 

Bye = (uh ap m4) One — 2m ge + m? i ee Q,(m) +- rae (m2 ie 4 Y,.)| 


Ge aan LORS 8 TP? a 


2. ay 
pe 25 Sal) — “gaat OP Jot Fa) Paya Ut tml) Fabian Day 


a 16 m? Ay tanh 1/2. Ae ays 2a 16 m? uw, coth 1/2 uy 
n Qs — mA (a4 — mi)? yn ak — mi (ud — mi)? ; 
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and ~ 2S — 22 
la = 76 — mi?’ Sn = ua — me’ 
1/2 
; P,(m) = | cosh m x sech > C,(x) dx , 
-1/2 
1/2 
Q,(m) = i sinh m x cosech =~ S,(%) dx , 
. —1/2 
1/2 
U,(m) = | (1 — 4?) cosh m x sech * C,(x) dx , 
= 1/2 
1/2 
Vi(m) = / (1 — 4 x?) sinh m x cosech > S,(x) dx . 
—1/2 


The quantities X,,, Ynq Zo. --- Zs, Jo--- Js which occur in the above expressions 
are constants whose numerical value is given in [6]; 6, is the Kronecker delta. 

For (33) to have a non-trivial solution, the determinant of the system must 
vanish; this provides a relation between 7, m and T, from which information 
about the onset of instability can be obtained. We shall approximate this 
infinite determinant by one of second order, which corresponds to expressing w 
as a linear combination of C,(x) and S,(x) only. The corresponding determinan- 


tal equation is then eB. 


4 Aes 
ek (34) 


By specifying 7, it is possible to find from (34) the relationship between T and 
m, and hence that between R(d/R,)? and m; moreover, we can find the 
minimum value of R(d/R,)"? (= Rerit(4/ R,)/?) and the corresponding m (= Merit): 
These results are shown in Table 1 and Figure 1. Table 1 also contains values 


Table 1 
Critical values of R(d/R,)1”” 
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Figure 1 
Curves of neutral stability. 
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a 
Figure 2 


The variation of [R(d/R,)"/?] 44 with T). 
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of B, (choosing A, = 1) at the critical points, and from these the streamlines 
and velocity profile at the onset of instability can be determined. 
Figure 2 illustrates the variation of R,,;,(¢/R,)/? with T, for T, > 0. 

_ The results obtained show that when 7; increases there is a decrease in both 
R.i,(4/R,)? and m,,;,. Results have not been obtained near 7, = 1707°8, at 
which value instability should set in with zero basic fluid velocity. This is 
because when V,, = 0 the stability equations are of such a nature that effectively 
yw is then being approximated by only one function. However, for the values 
of T, given in Table 1 it appears that the present approximation is quite 
adequate. The work might be extended to consider the effect of slight curvature 
of the walls on the stability of flow in a parallel plane channel under the 
influence of a magnetic field parallel to the planes. 
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Résumé 


On discute l’influence d’un champ magnétique circulaire sur la stabilité d’un 
fluide visqueux et conducteur de l’électricité en mouvement entre deux cylindres 
circulaires coaxiaux. Supposant que l’espacement entre les cylindres est petit par 
rapport a leurs rayons, on donne des résultats sur deux sortes de mouvement du 
fluide basique, produites par la rotation des cylindres et par un gradient de pression 
autour des cylindres. On démontre que le champ magnétique peut avancer ou 
retarder l’accrochage de l’instabilité, suivant les circonstances. 
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Uber die Wirbelstromdampfung 
bei ferromagnetischen Schichten 


Von Martin C, GuTZwILLeR, Ziirich?) 


Oe a ie tw nce am 


~- 


In der nachfolgenden Rechnung soll die Entstehung und die Wirkung von | 
Wirbelstrémen in einer einfachen Leiteranordnung untersucht werden. Die~ 
Situation kann in einem kartesischen Koordinatensystem (x, y, z) wie folgt 
beschrieben werden : 

In dem Zylinder ; 

Pty SR, |z|s2 (1) 
befindet sich ferromagnetisches Material, dessen Magnetisierung (M,, M,, M,) 
als homogen angenommen ist und in bekannter Weise mit der Zeit ¢> 0 
variiert, wahrend fiir ¢ << 0 die Magnetisierung nicht variiert. Uber und unter 
diesem ferromagnetischen Kérper befinden sich unendlich ausgedehnte leitende’ 
Platten, deren Position gegeben ist durch 


Sas |z| sd (2) 


und deren spezifischer Widerstand @ ist. Es seien die in den leitenden Platten 
induzierten Wirbelstr6me und die von diesen Wirbelstrémen erzeugten magne- 
tischen Felder innerhalb des ferromagnetischen Zylinders gesucht. 

Die Methode der Lésung lasst sich sofort auf Anordnungen der Leiter an- 
wenden, welche in komplizierterer Weise von der z-Koordinate abhangen. Dies 
wiirde zu verwickelteren algebraischen Ausdriicken fiihren. Aber die Ergebnisse 
waren im wesentlichen dieselben, insbesondere in dem wichtigen Fall, wo die 
vertikalen Distanzen, verglichen mit dem Radius R, klein sind. 

Die Maxwellschen Gleichungen fiir das magnetische Feld (H,, H,, H,) in 
den verschiedenen Gebieten lauten folgendermassen: 


ret; = Q, 


3 fine sliadpettee, (3) 
div(H +42 M)=0, 


wobei der Magnetisierungsvektor M als bekannte Funktion des Ortes und der 


1) IBM-Forschungslaboratorium, Adliswil-Ziirich. 
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Zeit vorausgesetzt wird; 


TOt ed —" Ol D 
fare ee eel eee arid: eel 0% 4 
fae elie aie j= | 4) 


das heisst, das magnetische Feld H ausserhalb der Leiter gehorcht den Glei- 
chungen der Magnetostatik; die Gleichungen fiir die Stromdichte (G,, G,, G.) 
in den leitenden Platten sind 

fo i Ce divG=0, 


ee fi ae es. (5) 


Ot Ce aa A div H=0, 


Die Randbedingungen an Unstetigkeitsflachen verlangen: 


Stetigkeit der Normalkomponente von H + 42 M, 
Stetigkeit der Tangentialkomponente von H, | (6) 


Verschwinden der Normalkomponente von G an der Oberflache 
eines Leiters. 


Um diese verschiedenen Gleichungen zu lésen, wird eine Laplace-Transfor- 
mation der zeitlichen Variablen und eine doppelte Fourier-Transformation der 
raumlichen Variablen (x, y) vorgenommen, das heisst der Operator 


co +co +00 


sere: a | i once dx dy (7) 


0 —0co -—0Oo 


wird angewandt. Es wird angenommen, dass die Integrale existieren und dass 
die transformierten Funktionen in der x- und y-Richtung geniigend stetig sind, 
um die Differentiationen nach x und y einfach als Multiplikation der Transfor- 
‘mierten mit 7 « und 7 erscheinen zu lassen. Da uns das fiir t < 0 bestehende 
magnetostatische Feld nicht weiter interessiert, so soll es sofort von H abge- 
zogen werden. Die zeitlichen Ableitungen von H erscheinen in der Transfor- 
mierten von H wiederum einfach als Multiplikation mit p. Sinngemass ist dann 
unter M nur der zeitlich veranderliche Teil zu verstehen. 

Es seien nun die transformierten Felder durch die entsprechenden kleinen 
Buchstaben m, h und g geschrieben. Ausserdem werde der symbolische Vektor 
f= (ia,7 8, 0/dz) eingefithrt. Dann werden die Feldgleichungen 


fx h}=0, (h+4am)=0, fir |z| <4: (8) 


pxbt=0, (fbl=0, fr > <l|z|<e, |al>>: (0) 
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[fxh]=-*g, (f,h)=0 
fir eae (10) 


elfxgl=—-fh, (f,8)=9, 


Die Randbedingungen an Unstetigkeitsflachen in der z-Richtung bleiben wie 
bisher. m ist als eine bekannte Vektorfunktion der Variablen «, 6 und # aufzu-_ 
fassen. Fiir grosse Werte von | z| muss das Vektorfeld h bei konstanten «, 6 | 
und # verschwinden. 

Die Beschrankung von M auf den Zylinder a + y? < R* soll erst spater in. 
der Rechnung beriicksichtigt werden. 

Die Gleichungssysteme (8), (9) und (10) kénnen immer auf ordentliche ee | 
ferentialgleichungen zweiter Ordnung zuriickgefiihrt werden, von denen nur 
die Gleichungen (8) nicht homogen sind. Mit den Abkiirzungen | 


y =o? + B, b= \/o8 + p+ ae ~ |/r+ +2? (11). 


ergibt sich in den verschiedenen Bereichen: 


2 
| 


Gg, 


(h,, hz) = an * (a, B) [1 sinhy z+ %, coshy z 


. " D 
es (am, + B m,)| , ¢ far |zi< 2? (12) 
h, = *,coshy z+ x,sinhyz—4am,, 
(h,, hy) = mi * (a, B) [A, sinhy z + A, coshy z] , | D 
fir = <|z|< 4; (13) 
h,=4,coshy z+ A,sinhyz, | 


(h,, hy) = : (x, B) . [uy sinhd z+ fy coshd z} 


| (14) 
h, = cosh 62+ we.sinhdz, 
. a a Eich nh 
(6a &) = 2 (= B, +0) PCy coshd 2 + pagsinh dz), rt 
ga), 
(h,, hy) = (x, B) [», siny z+, cosk yz], 


Ne Deacol 2 Al: (16) 
h, =, coshy z+», sinhy z, | 


Von den Randbedingungen ist eine bereits erfiillt worden, namlich g, = 0 
fiir | z| = a und fiir | z| = 0. Es ist daher aus (15) zu sehen, dass die Rand- 


Vol. XII, 1961 Uber die Wirbelstromdampfung bei ferromagnetischen Schichten Silbl 


bedingung fiir g hinreicht, um die Stromverteilung in den leitenden Platten 
ben zu machen, das heisst die Komponente g, zum Verschwinden zu bringen. 

Die in den Gleichungen (13) bis (16) auftretenden Koeffizienten A, u und » 
sind nicht gleich fiir z > 0 und z < 0. Darum werden sie im folgenden durch 
einen oberen Index + oder — gekennzeichnet, falls dies notwendig ist. 

Durch die Randbedingungen (6) werden nun zwischen den Koeffizienten 
x, A, und y lineare Beziehungen hergestellt. Das Verschwinden von h fiir 
grosse | z| erfordert, dass 


i, +e = 0, wm—m =. (17) 


Daraus folgt nach einiger Rechnung 


(a)xu, = 42m, ((a) cost > — (b) sinh a ) ; 


2 
47% D D oe 
4 é 
(by ag = — , (x m, + B m,) (= (a) sinh — + (b) cosh =) : 
wobei die Ausdriicke (a) und (b) gegeben sind durch 
fa) = e”* coshod (6 — a) + (“, coshy a + = sinhy a) sinhdé (b — a) , | 
(19) 
(b) = €”* coshd (b — a) + (~ sinhy a+ 5 coshy a) sinhd (b — a) . | 


Fiir die weiteren Rechnungen ist es eine grosse Hilfe, zwei spezielle Annah- 
men zu treffen, welche physikalisch sinnvoll sind. Erstens soll die Dicke D als 
klein gegeniiber den iibrigen Langen, die in dem Problem auftreten, angenom- 
men werden. Es sollen dann nur die niedrigsten Glieder in einer Entwicklung 
nach D beriicksichtigt werden. Zweitens soll angenommen werden, dass m, = 0 
ist. Dies kann dadurch gerechtfertigt werden, dass in einer diinnen ferromagne- 
tischen Schicht das entmagnetisierende Feld die Magnetisierung immer in die 
Ebene der Schicht zwingen wird. 

Es werden jetzt zwei spezielle Gréssen untersucht. Einerseits interessiert 
das in der Ebene z = 0 durch die benachbarten Leiter hervorgerufene magne- 
tische Feld h,,. Es wird aus dem Gesamtfeld h erhalten, indem man das Feld h, 
abzieht, welches auch in Abwesenheit der Leiter, das heisst im Falle 6 =a, 
besteht. Fiir dieses Wirbelfeld h,, erhalt man aus (12) und (18) den Ausdruck 


h = — 2% (A) (em, + Bm,) 2 hy), (20) 


a y? c2 


(h 


wx? 


wobei /i(y) gegeben ist durch 
1 e ’* sinh 6 (b — a) (21) 


My) = “5 78 cosh b (b — a) + [(d/y) sinhy a + (y/9) coshy a] sinh 6 (6 ea 
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. > 
Andererseits ist es wichtig, die Stromverteilung in den Leitern zu kennen. 


Man erhilt fiir die Stromdichte an der Oberflache | z | = a 
2 D 7 
(Ge &) = sat (— Bs +a) (am, + Bm) Esty) (22), 


wobei g(y) gegeben ist durch 


cosh 6 (b — a) + (y/6) sinh 6 (b — a) (23) : 
e’* cosho (b — a) + [(6/y) sinhy a + (y/6) coshy a] sinh 6 (b — a) ; | 


ely) = 


Aus (20) und (22) lassen sich die «Ausbreitungsfunktionen » fiir das Wirbel- | 
feld h,, in der Ebene z = 0 und fiir die Stromdichte g in der Ebene z = a be~ 
rechnen. Darunter werden in diesem Problem bei konstantem Parameter p 
diejenigen Matrizen mit zwei Reihen und zwei Spalten verstanden, deren Ele- 
mente von den Koordinatendifferenzen (x — &, y — 7) zwischen dem Beob- 
achtungspunkt (x, y) und dem Quellenpunkt (&, 7) abhangen und welche die 
Anregung, das heisst h,, oder g im Punkt (x, y) als Funktion des magnetischen 
Momentes m im Punkt (§, 7) angeben. 

Fiir das Wirbelfeld ist diese Ausbreitungsfunktion 


ea OR [fae apeinr By Oe Ay (24) 


Fir die Stromdichte g in der Ebene z = a gilt ein ahnlicher Ausdruck. 
In den Polarkoordinaten 


x—&=rcosp, y—yHn=rsingy (25) 


lasst sich (24) einmal integrieren, und es ergibt sich 


le,2) 


Tt 1 0 in 
~27bP fra tion (, §)- non (Soe Satin. 


Dabei ist eine bekannte Integraldarstellung der Bessel-Funktionen angewandt 
worden [1]?). 


Als nachster Schritt wird die Quellverteilung der Magnetisierung, wie sie 
durch (1) gegeben ist, beriicksichtigt werden. In den Polarkoordinaten 


&F=ocosy, n=asiny (27) 


wird der Ausdruck (26) iiber y von 0 bis 2 und iiber o von 0 bis R integriert. 
Wenn man ausserdem noch setzt 


%—=scosy, y=ssiny, (28) 


a Se 3 ue , ; ; 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 515. 
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so erhalt man weiter 


2iy a Gen 


Seg ee 


ry =|/s?—2s0 cos (p— x) +0, e2'¥ —e (29) 


Diese Formeln kann man in (24) einsetzen und das Additionstheorem von 
GraF [2] anwenden. Dann braucht man die bekannte Formel 


R 


R 
[odo Joy 0) = = hy o)[= RAE, (30) 


0 


welche aus den Rekursionsbeziehungen zwischen den Bessel-Funktionen folgt 
[1], und erhalt fiir die Laplace-Transformierte des Wirbelfeldes in der Ebene 
3 = 0 


aera meee 
wy 0 
(31) 


aos Z in2 VM 
x pi 5) | — flys) (Ge Se 4) hl) ea 
Ler sin2 yp —cosZp | M, 


Der Querstrich titber H und M soll andeuten, dass es sich um die Laplace- 
Transformierte beziiglich der Zeit handelt. In derselben Weise ergibt sich fiir 
die Stromdichte in der Ebene z= 4 


(32) 


= Y 
Oo -1 —sin2y cos2yp M 
x je s) ( — Jaly s) | | g() Bi 
+1 0 cos2y sin2y M, 


Der Parameter p der Laplace-Transformation steckt in M, und M , und 
ist ausserdem noch in A(y) und g(y) vorhanden, da letztere von 6 abhangen, das 


durch (11) gegeben ist. 
Ein Spezialfall lasst sich einfach behandeln, namlich 


a<R und 0<R, (33) 


wobei die Bedingung } < R dahinfallt, falls bereits in (21) und (23) die Lange 
b — a gross gegeniiber allen andern Langen angenommen wurde, was eine ent- 
sprechende Vereinfachung von (21) und (23) ergabe. Solange s mit R vergleich- 
bar ist, kommt die hauptsachliche Variation der Integranden in (31) und (32) 
von den Bessel-Funktionen her, aber nicht von h(y) und g(y). Diese letzteren 
Funktionen kénnen daher annaherungsweise durch ihre Werte fiir y=0 
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zuriickgefiihrt [3]: 
‘sg 1 (#\ tar A 
pate by ms ; 
A 
[a Daly») Id) =} (a) 
0 0 fir wu > A, 


und man erhalt aus (31) und (32) 


1 0\ /M 
Pe a) iy firs < R, 
H ee o 1/\M, 


S in 2 é 
ee Ae 40% p D 1(0) co. 2y sin4 py =. ar hie ee R, 
ee M,] \s 


sin2y —cos2y 


f 0 —1\ /m 
oe g(0) os furs <a ee : 
G ee +1 0/ \M, | 


7 
4 
Mi 
: 
| 
| 
; 


ee (36) 
—sin2y cos2 M ‘ 
: aes g(0) r a Me (=) furs Pik. 
e° cos2y sin2y/ \M,/ ‘° 
Die Grenzwerte von h(y) und g(y) fiir y > 0 sind nach (21) und (23) 
a snhd(b—a) | 
MO) eee nd i a) ad ow eed ee an 
Oyae cosh 6 (b — a) 
SiO) = Scand hs a) Oa ee Sal 
wobei jetzt 6 einfach gegeben ist durch 
_/4ap 
ja tS (38) 


Der Grenzfall 6 — a > oo wird auch noch aus (37) richtig erhalten, denn 
physikalisch bedeutet dies in vielen Fallen nur 6-! < b—a < R. Es ist nicht 
zu erwarten, dass die Formeln (35) und (36) im Bereich der Unstetigkeit, das 
heisst fiir Werte von s, wo|s — R| = a, noch irgendwelche Giiltigkeit haben, 
weil die Integrale (31) und (32) in Wirklichkeit tiber viele Schwankungen der 
Bessel-Funktionen mitteln und deshalb in der Nahe der Unstetigkeit sehr 
empfindlich gegeniiber langsam variierenden Faktoren sein miissen. Es kann 
deshalb vorkommen, dass die Vektorfelder nach den Formeln (35) und (36) bei 
s = R unstetig sind, obwohl dies nur bedeutet, dass sie im Bereich | s — R| = a 
sehr stark, aber stetig, variieren. Es wird sich allerdings zeigen, dass die durch 
(35) und (36) definierten Felder bei s = R in verstandlicher Weise unstetig 
sind, obwohl sie natiirlich dort keine stetigen Ableitungen haben. Es ist daher 
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verniinftig, anzunehmen, dass (35) und (36) im wesentlichen richtig sind. Zur 
naheren Erlauterung von (35) und (36) lasst sich zuerst sagen, dass das Wirbel- 
feld innerhalb und die Stromdichte oberhalb der ferromagnetischen Schicht, 
das heisst fiir s < R, homogen sind, und zwar das Wirbelfeld entgegengesetzt 
der Magnetisierungsanderung und die Stromdichte im rechten Winkel dazu. 
Ausserhalb der Schicht, das heisst fiir s > R, verhalt sich der Betrag der 
Felder H und G wie (R/s)? und schliesst stetig an die Betrage fiir s << R an. 
Im Bereich s > R erscheint die Schicht also wie ein magnetischer Dipol in 
zwei Dimensionen, was natiirlich eine Folge der Begrenzung der Wirbelstréme 
entweder durch die Dicke der leitenden Platten oder durch den Skineffekt 
darstellt. Die Richtungen der Felder (35) und (36) fiir s > R ergeben sich einfach 
aus folgender geometrischer Konstruktion: In jedem Punkt wird der Vektor 
(M,, M,) bzw. (M,, —M,) am Kreis s = konst gespiegelt. Es lasst sich dann 
kontrollieren, dass sich der Fluss der Feldlinien in beiden Fallen schliesst. 

Von besonderem Interesse sind die Grenzfalle einer diinnen leitenden 
Schicht, 6 > a, und einer unendlich dicken leitenden Platte, b > oo. Es muss 
dann in (35) und (36) nach (37) eingesetzt werden 


h(0)z (b—a), g(0)21 fir ba, (39) 
1 1 1 4 


wobei natiirlich 6 immer noch durch (38) gegeben ist. Im Falle (39) lasst sich 
die Laplace-Transformation der zeitlichen Variation ohne weiteres riickgangig 
machen, denn der Vektor (M,, M,) ist einfach mit p multipliziert worden, 
das heisst der Vektor (M,, M,) ist nach der Zeit differenziert worden. Im Falle 
(40) kann man 6 als die reziproke Eindringtiefe auffassen und fiir den héchsten 


Wert von # ausrechnen, fiir welchen (M,, M,) noch nicht vernachlassigt werden 
diirfen. Wenn sich dann ergibt, dass 6a < 1, kann man h(0) einfach durch 1/6 
-ersetzen und die Umkehrung der Laplace-Transformation ist dann in vielen 
Fallen méglich. 

Im allgemeinen Fall (37) ist zu bemerken, dass h(0) und g(0), ebenso wie (21) 


und (23), meromorphe Funktionen von p sind, so dass sich die Umkehrung der 
Laplace-Transformation als eine Summe tiber gediampfte Schwingungen dar- 


stellen lasst, falls auch (M,, M ,) meromorph ist. 
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Summary 


A ferromagnetic disc is placed between two conducting plates of infinite extent 
so that the axis of the disc is perpendicular to the plates. The magnetization in th ‘ 
disc is homogeneous and varies with time ¢ > 0 in a known but ‘held which tel 
The eddy current distribution in the plates and the magnetic field which the 5 | 
generate in the disc are represented with the help of a Laplace transform with” 
respect to time and a Hankel transform with respect to the distance from the 
center of the disc. In the special case where the distances perpendicular to the 
conducting plates are small compared to the radius of the disc, these two transforms | 
can be inverted. The eddy currents and the generated magnetic fields are then 
easily visualized. : 


(Eingegangen: 18. Februar 1961.) : 


Zut Kontinuum-Magneto-Gasdynamik 
Stationare rotationssymmetrische Stromung 
eines vollkommen leitenden Plasmas 


Adolf Busemann zum 60. Geburtstag 


Von Ernst HOLDER, Mainz 


Die klassische Gasdynamik [1]1) ebener oder rotationssymmetrischer Stré- 
mungen hat A. BUSEMANN [2] — fiir den Ingenieur wie fiir den Mathematiker 
gleichermassen lehrreich — dargestellt ; ich erinnere mich dankbar seiner Braun- 
schweiger Vorlesung. 

Heute, nach zwanzig Jahren, ist die Magneto-Gasdynamik [3] aktuell; ledig- 
lich deren Kontinuumstheorie haben wir im Auge. Ich méchte die Aufmerk- 
samkeit darauf lenken, dass auch hier die stationdren rotationssymmetrischen 
Strimungen eines vollkommen leitenden Plasmas in ahnlichen Schritten wie in 
der Gasdynamik behandelt werden kénnen: Ausgangspunkt ist das Mollier- 
Diagramm (A) fiir adiabatische Zustandsainderungen. Nachdem das Induk- 
tionsgesetz (Ff) die magnetische Induktion B aus Massenstromdichte ov und 
Massendichte @ zu berechnen lehrt, gibt von den Eulerschen Impulsgleichungen 
die auf das Azimut beziigliche Gleichung (£,) langs einer Stromlinie die Kon- 
stanz einer Linearkombination aus Drehmoment Q und azimutaler Kompo- 
nente B, der Induktion, und damit 2 als gebrochene lineare Funktion von 0. 
Sodann folgert man aus den beiden anderen Impulsgleichungen zuerst eine 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, S. 525. 
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durch B, erganzte Bernoullische Gleichung (B), weiter den Drehungssatz (D), 
die Erhaltung einer Linearkombination der azimutalen Komponenten € von 
Wirbelstirke und elektrischer Stromdichte k; das ist fiir die der Kontinuitat 
(C) geniigende Stromfunktion y die gewiinschte partielle Differentialgleichung 
zweiter Ordnung, eine Verallgemeinerung der Differentialgleichung fiir die 
Zirkulation einer Atmosphare [5b]. 

Nur im axialsymmetrischen Strémungsfall, wo auch die magnetischen 
Kraftlinien in den Ebenen durch die Achse verlaufend angenommen werden, 
scheint nach dem Vorbild des Busemannschen Druckberges, [2], S. 439, ein 
Druck-Magnetberg mit seinen Doppeltangenten die Sprungeigenschaften am 
Verdichtungsstoss darzustellen und die Figuratrix fiir ein Variationsproblem 
abzugeben [5a]. 


1. Strémungen mit adiabatischen Zustandsanderungen 


Das perfekte Plasma (das heisst elektrisch vollkommen leitend, reibungstrei, 
und nicht warmeleitend) mit der Zustandsgleichung zwischen Druck p und 
Dichte 9 = p; in Form des Mollier-Diagramms 


(A) p= pti s), 


Entropie s = s° = konst bei der vorausgesetzten adiabatischen Bewegung, En- 
thalpie i, ebenfalls pro Masseneinheit, kann sich in einer stationdren rotations- 
symmetrischen Bewegung befinden. Selbst dabei kann der Fermische «Accele- 
rationsmechanismus» von Ort zu Ort wirksam sein, der namentlich bei instatio- 
nairen Bewegungen interessant und technisch bedeutsam ist. 

Die zugrunde gelegten Differentialgleichungen [3] driicken neben der Diver- 
genzlosigkeit der magnetischen Induktion B und der Massenstromdichte 9 v 
das Induktionsgesetz und die Eulerschen Impulsgleichungen aus. 


tr), (Ff): V-B=0, Vxi[vx B)=90, 
(C), (E): V+ (ev) =9, wx vjxv=—V (E+ —— P+ 5, 


x VV x B)x B. 


Wir legen Zylinderkoordinaten x, y, Y zugrunde, in denen das Quadrat des 


Linienelementes 


dit dy+ydge?, g=y 


ist. SAmtliche Zustandsgréssen hangen nicht vom Azimut @ ab. y ist der 
Abstand von der Rotations-x-Achse. 
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Die Differentialgleichungen sind, in Komponenten nach dem orthonormier- 
ten Dreibein in Richtung der Zylinderkoordinaten”*) : 


2(nw—Cv) =—2. (T+ M) +7 GC kB), .. 4 


2(Cu-Ew)=- 5 (+H) a Et kA iC), (1.2) 


Lae : ' 
oP oral sol a : (1.3)% 
; : 1 
U(x, y) Kraftefunktion pro Masseneinheit, T= > |v|*, u (konstante) magne- 
tische Permeabilitat. Drehung V x v = 2 (§, 9, 5, elektrische Stromdichte 


il ena 
fa Ve Beg 
eae: (i, 7, R) 


| 
aa), 


2) In der Formel bei LizpMaNN-CoLE [3], S. 188, steht links in Cartes a 
Koordinaten 


CZ —a) = — —— 15) 
Al = aaa (1.5). 
in meiner Beziehung, wo die vektorielle Differentialform 
M =—BB+ 5 (B-B) do (1.6). 
und das dussere Differential 
dM = —dB (B- do) + (B- dB) do | 
= e ie Lae 
= —[do x dB) xX B= —[V x B] x B (dt) | ale 


sind. Dabei ist die Differentialform der Induktion (gleich in den Zylinder-Koor- 
dinaten) 


B= By, dy dy + By, dy dx + By, dx dy 


1.8 
= B' do, + B® do, + B* do, = A do, + B do, + C do = (B- do) | So 


mit 
B= Se B= 0 B= 51 a ; . Be (Bl By ce We Bic) | 
Gay dy) de “do, = ydpdx, do, aydudy=y do, i) 
do = (do,, do,, do) Oberflachenelement; (dt) = y dx ay dp Volumenelement. | 
Auf der rechten Seite kommt also mit (1.4) der Ausdruck 

J 1 aM 1 

fib Rea e these Sir [VX Bl xB: (1.10) 


Hier steht die «legitim-vektorielle Form» im Sinne von SOMMERFELD [4], S. 81, 326, 
die man ohne weiteres auf Zylinderkoordinaten umschreiben kann. 
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2. Magnetische Induktion 


Aus dem Induktionsgesetz (F) folgt mit den Komponenten wu, v, w der Ge- 
schwindigkeit v nach dem den Zylinderkoordinaten angepassten Orthosystem 
von Einheitsvektoren und mit den entsprechenden Komponenten der magne- 
tischen Induktion B = (A, B, C) fiir die im bewegten vollkommen leitenden 
Plasma allein vorhandenen elektrischen ‘Feldstarken 


Be—Cuo=-N., C#—AwsN,, Ava=Bbu=0 (2.1) 


% y 


mit einer Funktion N, fiir die 
uN,+vN,=0 (2.2) 


gilt, die sich also langs einer Stromlinie nicht andert. 
Die Kontinuitatsgleichung 


(C) (you), + (yor), =0 


kann man im rotationssymmetrischen Fall durch eine Stromfunktion y(x, y) 
mit der Eigenschaft 
yoeSy, =F, VOCs YK —P (2.3) 


befriedigen. Dann bedeutet (2.2), das heisst y, NV, — x N, = 0, dass N = N(y) 
ist. 
Aus der Divergenzlosigkeit der magnetischen Induktion (Fo), 


(Bos) x a (Bsi)y =0, 
die wegen (1.9) auch 


geschrieben werden kann, folgt ebenso 


yA=yBi=M,, y B= y B= —M,, (2:5) 


| y 
mit M = M(y). 
Also gilt mit M’(y)=m, N’(y)=m trae b= 5, b= ae = B,, 


C = B,/y = y B® 


(F) A=mou, B=mov, C=mow+noy. (2.6) 


3. Konstante Linearkombination von Drehmoment und azimutaler 
Komponente der Induktion 


Mit dem Rotationsmoment 
Q=yw= Ws (3.1) 
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sind im Impulssatz die Komponenten der Drehung V x v 
_— Laks Yaseen (3.2) 
Z é om y- 5 2 7] y , 6 x y 


und analog mit B, = y C die Komponenten der elektrischen Stromdichte 


1 
=—VxB, 
a 
ee oP ee en ee : 
ey Oy yy” nf wy Ox ay ) 
il 
k== (B,—24 
+ (B,- A, 
: 
Die 3. Impulsgleichung (1.3) hat die Gestalt : 
1 1 iS | 
se dy Oe $2, ti) ee (Ba), a eae oa 
integriert 
B 
(E;) Q— m 2s — Wy). 
Also gilt fiir das Rotationsmoment 2 die Beziehung 
m* m ee 
Q (1-7) -— ney, (3.5) 


die 


On Bit mney 


= WD ZWwe 
fu —m? 9 


(3.6) 


durch 0, y, y ausgedriickt gibt. 


4. Durch Induktionskomponente ergianzte Bernoullische Gleichung 


Nun kann noch die Dichte 0 = p;, namlich i, vermége einer Bernoullischen 
Gleichung (B) langs jeder Stromlinie durch die Komponenten #, g der Strom- 
dichte, das heisst die 1. Ableitungen von yw ausgedriickt werden. Dazu benutzen 
wir die kovariante Webersche Form der Bewegungsgleichungen (1.1), (1.2) mit 


den Drehungskomponenten (3.2) und erhalten durch Faltung mit w, v mit der 
Bezeichnung 
0 0 


—D,(t+%4+2 tu) + Sv: [jx Bl = (4.2) 
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Hier ist der Warmeinhalt i= u + p /o, u die innere Energie pro Masseneinheit, 
und fiir die adiabatische Bewegung der 1. Hauptsatz beriicksichtigt, mit der 
absoluten Temperatur T = —p,/e 


pD,—+D,u=TD,s=0. (4.3) 
Wegen 
1 ; ie 
<0 [jx Bl= 7 j-(B x0] 
1 > (4.4) 
en to Rate Amma DAB) 
R fallt heraus, kommt also integriert die Bernoullische Gleichung 
(B) T+UFE+ = = hy) |, 
beziehungsweise 
1 . Q ; nw = n ih 
4 t+) +i+2(se + F) er YN. (4.5) 
heraus. 


Aus der Bernoullischen Gleichung (8) folgt wegen (3.6) nach kurzer 


Rechnun 
= (v2 + v?) = ip —t— R (i) (4.6) 


mit 


2 
(mi+ ny) 2H —m Pi ee 
Ri, X,Y; 2) = 242 (u —m* p,)? Pj 
und aus (2.3) 
y2 p22 [ip -i- RY] =P +e, (4.8) 


wird dies nach i aufgelést: i= 3(p? + g?; x,y, 2), So ist (wv, v) durch (p, 9) 
ausgedriickt : he oes 

q : = re 4.9 
“= 5 >, OTR)? * y PB; (S(p? + 2") Ge 


5. Linearkombination der azimutalen Komponenten von 
Wirbelstarke und Stromdichte 
Endlich bekommt man wieder den Drehungssatz durch Anwendung der 
Differentialoperation 


0 0 a ee Sat 
Djeme e t bzw. ae a yo (uw? + v2)’ ( ) 
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Q 
(w+ v9) 26-20 (Eut nr) = (+r) 20——, D, 2 
) 1 
=D. (4 ae ry eevee ie. 
wo die 3. Komponente 
K,=(vx [jx B),=4[f x Blh—v[j x Bh 
=u(kA—iC)—v(jC—kB) 
= —©_ (— o(B,), + u(B,),) + & (w A +0 B) 


wy 
in 


aT 


So erhalt man den eee 


t= % [Dy (2 mt) + mes] ioemiece 


me) oD, Zt + mk (ut +P) @. 


e (5.4). 
re ae cg rs TD, { 
falls adiabatische Schichtung | 
pD,-- + D,u=TD,s (5.5) 

besteht. Oder mit 

Q= 1+ m By 
1 
20—mk _— 2 (y | m' Bs\ _ pi —_ Bs ree 

(P) ye = ( 53 be ( le )+Ts ’ bo) 


wahrend rechts nur Terme in den Ableitungen 1. Ordnung auftreten, steht 
links der negative Differentialausdruck 2. Ordnung in »p: 


Betty hice)” Yo m0 Oo wy  m (— Ta 0 yy ee) 
u 


Ov yo oy y Ox 
} é 4 ‘ (5.7) 
Place a ay ce 
0x L\e@ hing oy L\@ ew] y ue Var fe 


Der 1. Term ist 


r(0- M9) g1--([0- stale), 


2 


Vol. XII, 1961 Zur Kontinuum-Magneto-Gasdynamik 523 


Il 


3 2 
nit he Met LEE 59) 
= | By WV, LO ist die Alfvén-Geschwindigkeit fiir eine «inkompressible» 
Transversalwelle, die Drehflache normal zur Komponente By = (A, B) des 
Magnetfeldes in der Meridianebene ist. Wir erhalten im ganzen eine elliptische 
Differentialgleichung 2. Ordnung in yp. 


6. Variationsproblem fiir die axialsymmetrische Strémung 
eines Plasmas 


Im Fall einer axialsymmetrischen Strémung, wo auch die magnetischen 
Kraftlinien (identisch mit den Stromlinien) in der Ebene durch die Achse 


verlaufen, 
DQ=yw=0, B=yC=0, n=0, 1=0, (6.1) 


zeigen wir noch, dass die Differentialgleichung (D) die Lagrangesche Gleichung 
eines Variationsproblems ist. Wir nehmen als Grundfunktion 


f=t%y2PD=f+h, (6.2) 
A Nee AE Bel Sa pale (6.3) 

wo 
fo=y (pte (+r) (6.4) 


die alte Grundfunktion fiir das Variationsproblem der rotationssymmetrischen 
Bewegung eines Gases ist, aufgefasst als Funktion von x, y, z; p, 9, wie oben 
bei (4.9) allgemein ausgefiihrt ist: 


ha Sere + (eae (6.5) 


Bei einer ebenen Strémung fehlt hier und im folgenden der Faktor y. Dazu 
kommt der vom Magnetismus herriihrende Zusatzterm 
See ape gt 6.6 
f= +e). (6.6) 
Dann wird 
m2 


m2 
s 7, = 6.7 
FS=},=—-h, tae?) 7 “fate. 6.7) 
Bei (6.1) gilt die Bernoullische Gleichung (B) mit R(i) = 0, und (6.3) gibt 
den alten Zusammenhang zwischen (u, v) <> <= (f,9) — wie in [5a], S. 373. 
Wegen der Reziprozitat der Legendre-Transformation bekommt man 


—fop =, —fyg=—, (6.8) 


Auf der einen Seite ist die Divergenz 


af y Gfg «a dn 
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somit 

: fe 1 6.9) 

_ 


dx‘ dy dx dy 
oy) Ou m [ 0(m ev) O(meu)\ m’ 3 e 10) 
a ey ple Ox Oy i ra cea ( 
ste BP 5 gh 
2o+m ee ) : 
Auf der anderen Seite berechnet man aus (4.6) 
1 ‘4 
t= feet he=y (G2 +P. Ge) - Ge +) 
. (6.11) | 
= ye lh (y) — Ts'(y)] — = +). 
Das Variationsproblem 
| [te y, 2; p, q) dx dy — stationar (6.9) 
fiir die Stromfunktion 
ee _ oY _ oy 
z= y(x, y) , caer a e Bye 
hat die Lagrangesche Gleichung 
afy Ug _ 
paar a =f,. (6.12) 
Das ist 
=20 +m k= yo (ho =T s)i, (6.13) 


unser Wirbelsatz (D). 


7. Druck-Magnetberg als Figuratrix 


Gleichbedeutend mit (6.12) sind die Integralrelationen fiir jedes Teilgebiet 
T der Meridianebene und seinen Rand 0T 


[ody —nde= [G,dxdy, [Ede— Cav = [G, dx dy, 
oT iy 


(7.1) 
[ndx— 8 ay = [G. dx dy 
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mit dem Druck-Magnetberg als Figuratrix 


C= ey ty, 2) =) Pi, — U1, (7.2) 
vole bE se) 


: 2 : il 
wo in p und p, das Argument i = ty — > (u? + v?) steht und vermdge 


(u? + v?) [1 — a P; (i, os = (u2 + »))] = 4 72 


u2 + v2 auf €2 + 72 umzurechnen ist. 
Fiir einen schmalen Streifen zu beiden Seiten eines Verdichtungsstosses S 
folgen daraus die Sprungrelationen 


(E]: [ny] : (f] = dx : dy : de (7.3) 
fiir den Sprung 
(Ghent ae 


Diese besagt nach [5c], dass die Zustandsanderung E — £ bei einem Stoss durch 
eine Doppeltangente dargestellt ist, die von einem Pankt. (6°, 7°, ¢°) der zu 
s = s° gehorenden Figuratrix zu einem Punkt (é, n, €) einer zweiten Flache 


C= (6,;s) mit s>s? (7.4) 


fiihrt, in dem diese von der Tangente ebenfalls beriihrt wird, wie es BUSE- 
MANN [2], S. 439, angegeben hat. 
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Raumes). Arch. Math. 5, 510-521 (1954). 

[5c] E. HOLDER, Die Differentialgleichung fiir die Stromfunktion dey stationaren 
yotationssymmetrischen Zirkulation einer adiabatisch verdnderlichen veibungs- 
freien Atmosphiare tuber dey Evde, Meteorolog. Rundschau 7, 449-450 (1948). 
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Die in [5a] und [5b] dargestellten Variationseigenschaften der klassischen 
Gasdynamik und ihren Zusammenhang mit den Sprungrelationen am Verdich- 
tungsstoss habe ich zuerst in meiner unverdffentlichten Lilienthal-Preisschrift 
entwickelt, die 1943 in der Luftfahrtsforschungsanstalt Braunschweig entstan- 
den ist, im Aerodynamischen Institut von Professor BLENCK. 


Anmerkung bei der Korrektur: 


Erst nachtraglich habe ich (auf ebene Potentialstrémungen beziigliche) Va- 
riationsbetrachtungen bemerkt, bei A. BuSEMANN, Problems of Flow Existence 
and Stability in Magnetohydrodynamics in ‘The Magnetodynamics of Conducting 
Fluids’, ed. by D. Bershader, Stanford Univ. Press Calif. 1959, p. 3-16; ferner 
die auf einem isoperimetrischen Variationsprinzip von S. CHANDRASEKHAR be- 
ruhende Formulierung der Grundgleichungen gewisser rotationssymmetrischer 
Strémungen (unsere s’ = n’ = m’ = 0, h’ = —I’ n/m = konst) bei L. WoLTJER, On © 
Hydromagnetic Equilibrium, Proc. Nat. Acad. of Science 44, (1958), p. 833-841. | 


Summary 


In the case of rotational symmetry for the steady motion of a perfectly conduc- — 
ting plasma a differential equation of 2nd order for the stream function is given — 
by completing BERNOULLI’s equation, the equation for the moment of rotation and ~ 
the equation for the vorticity respectively by the azimutal component of the 
magnetic field and the azimutal component of electric current density. 


(Eingegangen: 14. Marz 1961.) 


Uber die Grenzwahrscheinlichkeiten einer speziellen 
homogenen Markoffschen Kette’) 


Von Hans-JoacHim TOpFEr, Berlin, Deutschland 2) 


Bei gewissen physikalischen Fragestellungen wird man auf Markoffsche Ket- 
ten mit endlich vielen, sagen wir N + 3, verschiedenen Zustanden By, «<> Enos 
gefiihrt, die durch die folgende Ubergangsmatrix beschrieben werden: 


1) Die Anregung zu dieser Fragestellung verdanke ich Herrn Dipl.-Phys. J. Biersack vom 
Sektor Kernchemie des Hahn-Meitner-Instituts fiir Kernforschung Berlin. Herr Brrrsacxk hat das 
in diesem Zusammenhang auftretende mathematische Problem in anderer Form behandelt. Unsere 
Zusammenarbeit hatte zum Ziel, eine fiir die numerische Bearbeitung durch ein elektronisches 
Rechengerat geeignete mathematische Form zu finden. 


*) Hahn-Meitner-Institut fiir Kernforschung, Sektor Mathematik, Berlin-Wannsee. 
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TeV OM as: O35 en 0 
ee a r 
Oe Pe aa 
P=||?,;||= : 5 as ; ig=1,...,N+3. Q 
DRO AG Os # 
Onipsonige 
Oey ol U 
einn gel see piesa" Ue. Us 1 
Man denke etwa an ein Elektron, das endlich vieler Energiestufen (Fy, ... , Ey +2) 


fahig sei und sich in einem «zufallsgesteuerten» elektrischen Feld befindet. Ist die 
Wahrscheinlichkeit dafiir, dass das Elektron zwischen zwei Beobachtungen auf die 
nachst-niedrigere, bzw. nachst-héhere Energiestufe gebracht wird, p, bzw. q, und 
wird es wahrend der gleichen Zeit mit der Wahrscheinlichkeit 7 durch Stossvorgange 
vernichtet (Zustand Ey +3), so gelangt man zu der angegebenen Ubergangsmatrix, 
wenn man annimmt, dass die Energiestufen FE, und Ey +2 sowie der Zustand Ey +3 
stabil seien, das heisst, einmal erreicht, nicht mehr verlassen werden kénnen. 


Die Elemente f;; der Matrix P geben die Wahrscheinlichkeit dafiir an, dass 
sich das betrachtete System nach einem «Versuch» im Zustand E;, befindet, 
unter der Bedingung, dass es vor dem Versuch den Zustand E; besass. 

Bekanntlich gilt fiir solche Ubergangsmatrizen notwendig 


Diy = 1, (2) 


wir werden also 
pagar= 1 (3) 
voraussetzen. 
Die Matrix der Wahrscheinlichkeiten fiir den Ubergang des Systems vom 
Zustand E; in den Zustand E;, nach einer Kette von » Versuchen ist nach MAR- 
xorr P”. Uns interessiert hier der Grenzfall einer unendlichen Kette von 


Versuchen. 
Die zugehérige Ubergangsmatrix bezeichnen wir, wenn sie existiert, mit 


p® =lim P*. (4) 
Die Elemente der Matrix P® = || p{°° || nennen wir Grenzwahrscheinlichkeiten. 


Mit Hilfe des folgenden Satzes aus der Theorie der stochastischen Matri- 
zen) lassen sich diese Grenzwahrscheinlichkeiten berechnen: 


8) Eine quadratische Matrix heisst stochastisch, wenn ihre Elemente nicht negativ und alle 
ihre Zeilensummen gleich eins sind. 
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Ist P die Ubergangsmatrix einer schwach regelmdssigen Markoffschen Kette, 
das heisst, besitzt P keine von 1 verschiedenen charakteristischen Wurzeln mit dem 
Betrage 1, so existiert P®. Sei y(A) das zu P gehorige Minimalpolynom und 


Cl) = (AE — P)* (0) 5) 

die reduzterte adjungierte Matrix zu P, so st | 
oo _ C(t) a. | 

IP = SS at 6 | 

v'(l) } oN 


ae 


(Fiir den Beweis vgl. [1]*) 
Um den Satz anwenden zu kénnen, miissen wir nachpriifen, ob die Voraus-_ 
setzungen in unserem Falle erfiillt sind, das heisst, ob P schwach regelmassig ist. — : 
Dies trifft in der Tat zu, denn wir k6nnen folgendes beweisen: | 
Das charakteristische Polynom der Matrix P ist | 


A(a) = (A 1)? Dy(Q) (7) 
worin Dy(A) durch die Rekursionsvorschrift 
D(A=1, D,A)=4, 


(8) 
Dy+2(4) =4*DayilA) —P9D,(A) n=0,1,.. 


bestimmt ist. Samtliche Wurzeln von Dy(A) sind einfach und reell und gegeben 
durch 


— k t 
A= 2) 04 Oe 4 ess NG (9) 
Beweis: Das charakteristische Polynom von P ist durch die Determinante 
A(d)=|AE-—P| (10) 
definiert, also 
Ae 0 oe Oger Oy a See ee 0 
—p hk -q O —r 
0 =-~f A —q 
A(A) = (11) 
—p A -q —»r 
0 0O A-1 0 
0 ee ee A ans aes Tl 


4) Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 535. 
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Hieraus sieht man zunachst, dass 


iy OEE Ok ei er ae 
—p A -gq O 
0 -p~ A -q 
Ma) =(4—1)3 | . Pleas (12) 
=p hg 0 
Oi =p Arg 
Orme 3 : : F . 0 0 -p A 


gilt. Die hier auftretende N -reihige Determinante nennen wir Dy/(A). Definieren 
wir 

D(A) =1 
und 


so besteht fiir beliebiges n = O die Beziehung 


Ae Gg, OO —$ —¢ Vir, 0 
=p 0) 

Di +2 (A) a 0 D ald) = A D+1(A) fe q f D,(A) 
0) 0) 


= 41D, 414) — bg D,() » 


das heisst, unsere Behauptungen (7) und (8) sind fiir beliebiges n = 0 richtig. 
Um nun (9) zu beweisen, stellen wir fest, dass zwischen den Polynomen 

D(a) und den Tschebyscheffschen Polynomen zweiter Art U,(x) mit dem 

Hauptkoeffizienten 1 ein einfacher Zusammenhang besteht. 


ZAMP XII/34 


, 
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Die Polynome U,(x) folgen der Rekursionsvorschrift (vgl. [2]): 


(13) 


Ersetzen wir in dieser Formel x durch A/2//pq und multiplizieren auf beiden 
Seiten mit (2 |/p q)"*2, so entsteht: 


(2) 64)" Case (aye) =A VPA) Oooa (Siew) 


~paevpay'0a(st=). ca 


- 
| 
Mithin ist die Rekursionsvorschrift (8) erfiillt, wenn wir 
| 


(yea) ales Osi 
setzen. 


Die U,(x) bilden ein Orthogonalsystem mit der Belegung V1 — +7 im Inter-— 
vall (— 1, 1), folglich bilden die Polynome D,(A) ein Orthogonalsystem mit der 
Belegung //4 pq — A? im Intervall (— 2 /p q, 2 /p q). Das bedeutet aber nach — 
allgemeinen Satzen iiber orthogonale Polynome, dass samtliche » Nullstellen- 
von D,(A) reell und einfach sind und im Intervall (— 2 |p q, 2 ’p q) liegen. 

Da man die Nullstellen der Tschebyscheffschen Polynome kennt, sie sind — 
namlich 


y= COS eet kR=172).32.05 (16) 
sind die Nullstellen der Dy(A) wegen des Zusammenhangs (15) leicht anzugeben 
und berechnen sich nach Formel (9). 

Weil nun, wegen r => 0, $+ q <1 ist, was bekanntlich auch 2 Veq S11 
nach sich zieht, ist also 1 die einzige Wurzel von A(A) vom Betrage 1, das heisst, 
P ist schwach regelmiassig. 

Wir schreiben noch das Minimalpolynom y(a) von P auf. Es ist: 


p(A) = (A — 1) Dy(A) - (17) 


Den Ausdruck (6) fiir P® wollen wir mit Hilfe der zu P adjungierten Matrix 
B(A) umformen. Wegen B(A) = (A E — P)-1 A(A), ist 


CA) = Fay 9) = aaa. (18) 
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Da 
i y'() = Dy(d) + (A 1) D4), (19) 
wird 
y'(1) = Dy(1) + 0. (20) 
Somit finden wir 
P® = lim, 24 (21) 


it (A— 1)® Dy(1) * 


Wir berechnen zunachst von P® ein allgemeines Element pir) Mit. cya, Ae 
N + 1. Dazu benétigen wir das entsprechende Element 0;, der Matrix BA), 
das bis aufs Vorzeichen mit derjenigen Determinante iibereinstimmt, die aus 
A(A) durch Streichung der i-ten Spalte und j-ten Zeile hervorgeht. Da diese 
Determinante aber, wie man sofort sieht, den Faktor (A — 1)8 enthalt, wird 
wegen (20) und (21) 


A= 0 fir 1,7=2,...,N+1. (22) 
Ist 


| 1 | 1 
i=IN4+2, aber 7#(N+2; 
N+3 Bee 


so entsteht aus A(A) durch Streichung der entsprechenden Spalte und Zeile 
eine Determinante, die eine Nullzeile enthalt. Folglich ist allgemeiner 


po) =O fir j=2,..,N+1. (23) 
Wir betrachten nun nacheinander die Fille j = 1,7 = N + 2 und 7 =N + 3. 
1. 7 = 1: 
Ist 7 = 1, so findet man das entsprechende Element von BA) zu 


by (A) = (A — 1)? Dy (A) , 


also ist 
pm) =1. (24) 


Ist dagegen i = N + 2 oder t= N + 3, so entstehen wieder Determinanten, 
die eine Zeile mit lauter Nullen besitzen, das heisst 


Pn ae = Neat =0. (25) 


Sei nun 1=2,...,N +1, dann erhalten wir aus A(A) durch Streichen der 
ersten Zeile und i-ten Spalte das entsprechende Element von B (A) zu: 
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—piaiHg0 ae 
0 -p A -q 
; 0 
oP Pad 
eine 
0) i Ora p="G 
b,=(—1i*A—-1) calle alll aaa FG 0 Oh? 
—pr-q 
0 Panag 
Or 2-5. Or =prk 
ee Se 
N—-1+4+1 
das heisst } 
by =P? (A= 1) Dy_ i) (27) 
woraus 
pe) = pit See) (28) 


folgt, oder mit Hilfe der Tschebyscheffschen Polynome ausgedriickt: 


5 : Un - ( ; WT ) 
q 
eles aoa (29) 
2Vpq 

209 = N+ 2: 
Wie oben findet man 

OE Ae = Py awee =0, (30) 

PN v2,N42 =1. (31) 


Sel nun t= 2,...,.N +1, so erhalten wir aus A(A) durch Streichen der 
N + 2-ten Zeile und 7-ten Spalte 
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Re = 0 ie og 0 
ap Ag 
is 0 
Pan 9 
(ost 0p A 
es 
; om -p-g 0 0 0 
me (1 N+2-4(7 __4)2 t—2 
+2 ( ) (A ) Z =O 0 ) 
=p = 
-pi -q 9 
Oia pt A Tg 
N—i14+2 
also: 
4 Fi bi noe = a (A— 1)? D,_ (A) - (33) 
Damit wir 
(0) _ ,N-i+2 De—a(1) 
Pinte = 47 “Dwi() (34) 


oder anders geschrieben: 


i 1 

a 1 ——\ N-i+2 Clee | 

(is > oat (Vz) ‘ ie (35) 
N —= 
( 2VPq 


a. gee St 


Wie oben findet man 


(00), (0) = 
1,N+3 PN+2,N+2 — 0, 


ie (37) 


(36) 


PN) s, N+8 — 
Der allgemeine Fall erledigt sich am schnellsten durch die Bemerkung, dass die 
Grenzmatrix stochastischer Matrizen wieder stochastisch ist, das heisst samt- 
liche Zeilensummen von P® miissen gleich eins sein. 
Daraus folgt: 


bee aaa ye 120 NY ey 


(32) 
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oder $),.,=1 pi-1 Dn—i+1(1) + gN-#+? Di-2(1) (39) 
1,N+3 Dwn(1) 


Eine besonders einfache Form nehmen die Grenzwahrscheinlichkeiten in dem 
Fall an, dass entweder # oder g verschwindet. 
Dann ist namlich auf Grund von (8) 


Di) =i fir pg=0, (40) 
so dass fiir = 0 
imeem gS.27 eal OM (41) 
fair ¢,= 0 : 
pS =p", PENse = (42) 
gilt. 


‘Andererseits wird der Fall haufig sein, dass p = g ist. Verschwindet zusatz- | 
lich y, so kommen wir wieder zu sehr einfachen Ausdriicken. Da dann p=q=1/2 { 
ist, bendtigen wir, wenn wir die Formeln (29) und (35) anwenden wollen, die 
Werte der Tschebyscheffschen Polynome an der Stelle x = 1. 

Man kann sehr leicht ausrechnen, dass ; 


Cte (43) | 


ist, womit wir aus (29) erhalten: : 


(eae ergy ae (44) 
und aus (35): 
co 1) — 1 : 
(ha aE (45) | 


Diese Formeln gestatten eine einfache und anschauliche Deutung: Wir 
denken uns, das System erstrebe den Zustand Ey, . Dann bendtigt es, um 
vom Zustand EF, aus auf geradem Wege sein Ziel zu erreichen, N + 1 Schritte. 
Befindet es sich bereits im Zustande E;, so hat es gewissermassen schon 7 — 1 
Schritte auf sein Ziel hin vollbracht. Die Wahrscheinlichkeit, das Ziel zu er- 
reichen, wird dann durch den Quotienten aus der Zahl der bereits zuriickgeleg- 
ten Schritte und der Zahl der benétigten Schritte gemessen. 

Wir wollen zum Abschluss noch eine fiir die numerische Auswertung 
zweckmiassigere Schreibweise der Formeln (29) und (35) angeben. 

Zu diesem Zwecke erinnern wir daran, dass die Tschebyscheffschen Poly- 
nome U,(x), die durch die Rekursionsvorschrift (13) offenbar fiir alle x erklart 


sind, im Innern ihres Belegungsintervalls (— 1, + 1) eine zweite Darstellung 
besitzen. Es ist namlich: 


= 1 sin n 1 arcco 
n(x) ae : = sacks 


gh Tas | wo lS hee (46) 


= 
if 
S 
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Die in dieser Darstellung zunachst geforderte Beschrankung auf den Bereich 


|x| = 1 kann nicht wesentlich sein, da U,(x) einPolynom ist. Die entsprechende 
Erweiterung fiir x > 1 ist leicht anzugeben zu: 


= ke = sinh (n As 4) Arch # 
Ve-1 

Da mit Ausnahme des schon behandelten Falles # = g = 1/2 stets 1/2 VPq =1 

ist, miissen wir (47) benutzen und erhalten: 


po =(/2) 


yea (47) 


: ; 1 
sinh oe —i-+ 2) Arch ae He 


? 


; if 
hie 1) Arch ——— 
sin { + 1) Arc Ata 


: 1 
ees pe sinh @@—1 Arch ———— 
mV el 
HN+2 p 1 
sinh {(N + 1) Arch _——| 
\ 2Vp4q J 
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Summary 


Simple expressions are given for the asymptotic probabilities of a special 
homogeneous Markov chain which is important in some physical problems. For 
this purpose a theorem about stochastic matrices is used which allows to compute 
the «co — th power» of such a matrice, if there exists no eigenvalue A with 4 + 1, 
but | 2| = 1. This condition is easy to verify, since a full survey of all eigenvalues 


can be obtained. 


_(Eingegangen: 27. Marz 1961.) 
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The Deflection and Buckling of a Beam-Column with a One- 


Sided Constraint ; 
By Epwarop E. Zajac, Murray Hill, New Jersey, USA}4) 


Introduction | 


We consider the planar problem shown in Figure 1: A hinged beam with a 
central vertical load Q is compressed by a load P. However, the beam is_ 
constrained, say by a floor, so that it cannot deflect in the direction of the Q. 


f 4 p 
Vs 
(Q-n)/2 N (Q-N)/2 
Figure 1 Figure 2 
Constrained beam-column. Deflected shape for P = Py. 


Without the load Q, one knows that deflection can occur if the load P 
exceeds the Euler buckling load P, for the pin-ended beam. But how much is 
this load raised by the presence of Q? One might attack this problem from a 
static viewpoint in the standard way by asking for the lowest load P at which 
there is possible a nontrivial equilibrium configuration infinitesimally close to 
the straight configuration. This load and configuration are immediately 
apparent; namely, that shown in Figure 2, with a load P, equal to the buckling 
load of a beam of half the span length hinged at one end and built-in at the 
other. The load P, is 8-18 times higher than P, regardless of how small Q is. 
Thus, using this approach, one would conclude that he could raise the buckling 
load by more than eightfold by applying an arbitrarily small transverse load Q. 

On the other hand, for very small Q one expects that deflection of the rod 
will be possible for a load P near the pin-ended Euler load P,. Hence, at the 
outset it appears that this problem may not be completely elementary. Indeed, 
we find first that the problem has several interesting features, which we 
investigate both by elementary beam theory and elastica theory. Secondly, 
we find that this problem is another example of what has been pointed out by 
ZIEGLER [1]*), namely that the classical elastic stability criteria may lead to 
difficulties. In this case, the difficulty is of a different sort than that discussed 


1) Bell Telephone Laboratories, Incorporated. 
*) Numbers in brackets refer to References, page 546. 
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by ZIEGLER. We find that a meaningful stability criterion is a function not 
only of P and Q but of the amount of the disturbance as well. 
Solution by Elementary Beam Theory 


By symmetry we have as the differential equation and boundary conditions 
for elementary beam theory 


ad = 
eG Le aeeleS (1a) 
y(0) aa =0, (1b) 
ae (10 


where JN, the reaction force of the floor at midspan (Figure 2), is zero for 

y(L) > 0. In addition we have the constraint condition (x) = 0 over the 

beam. The solution of Equation (la) for the boundary conditions (1b, Ic) is 
(Q —N) P 


2 pe eds a 
is Pa (q+ rae q rs (2a, b) 


The constraint condition y(x) = 0 gives 


_ sing* x 
i* cosq L 2 0. 


The only positive value of cos g L for which this is satisfied is cos g L = 1, or 
Sela A647, .<: , this corresponds to dy/dx = 0 at x = 0 and is of secondary 
interest. For negative cosqg L, y(x) = Oif (a/2) SgLs &, where & = 4-493 is 
the first root of € = tan. The value of g L = 2/2 corresponds to the Euler 
load P, (pin-ended beam), while g L = — corresponds to P, (pin-ended, built-in 
beam of half length). For gL < &, y(L) > OandN=0. AtgL =& (P= P,), 
y(L) = 0 and N may have any value between zero and Q. For (7/2) S<qL <a 
the deflection is of ‘first mode’ type, with an inflection point only at the hinge. 
For z < gL <&, it is of ‘second mode’ type, with an inflection point at the 
hinge and one between the hinge and the central load. 

For q L > é the situation is more complicated. First, there is the case where 
the beam touches the floor at two points between hinges, shown in Figure 3. 


Q 


|x| 


Yn, We WY YU Wye 
(Q-n)/2 ING N/2 (Q-N)/2 


Figure 3 
Possible deflected shape P > P3. 
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Taking the deflection, slope, and moment to be continuous at the touchdo 
points, we find 


MG—M) (gx — Si*), 08x81, 


id oe Zire cos & 
G=x2—F, 
Tye q(L—l Late 
y= gly [tn(tE=D) mat at sina], Ty 


L—lI 
(l-§, Q-N=—o Mle 008) 


The constraint condition y(x) = 0 is satisfied only if Q — N = 0, which pia 
tan [¢ (L — 1)/2] <0. The lowest values of ¢ L that satisfy this condition are- : 
given bya +é<qLl <2 +. Likewise, we may have higher modes of this: 
type with more than two touchdown points. | 
On the other hand with g L > &, the beam may have no touchdown points 
between hinges. This is possible if cosq L is negative and | cosq L | = | cosé |. 
The next values after gq L = & to satisfy this condition are in the range : 


4n—&SqlLs2n+. 


Counting the hinge these contain three or four inflection points in0 S * SL. 
Again higher modes of this type, with still more inflection points, are possible. 
Further, in the entire range 47 —€ Sq L S22 + € one may have either zero 
or two touchdown points between hinges. 

Hence elementary theory gives several interesting results. It indicates that 
deflections are possible for P = P, inarange of P less than the first ‘nontrivial 
load’ P,, but with the physically untenable result that y grows without 
bound, Equation (2a), for P +> P,. This is not surprising. For the first order 
elementary theory the load P, just equilibrates the bending moments at any 
cross section for any deflection if @ = 0. Hence, if Q + 0 and P is in the neigh- 
borhood of P,, we can expect that higher order effects will have to be considered. 
We note also that the spectrum of P’ s for which equilibrium states are possible 
is neither discrete, as in a linear buckling problem, nor continuous, as in 
ordinary beam deflection. Rather it is disjoint. Equilibrium configurations 
exist for (7/2) Sg L < & but donot exist for§ < gL < 2a, 0r2a<ql<a-+é, 
etc. Further, for (7/2) <q L < € it is seen that the deflection at midspan, 


Q 


So Page! ee (q L —tang L), 


surprisingly decreases as P increases. This may be explained as an arching 
effect; as P increases a shallower and shallower arch is required to maintain 
equilibrium. 
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Some of these effects are due to the omission of higher order terms in 
elementary beam theory and some to the constraint y(«) = 0. For a closer 
look, we therefore turn to elastica theory, where we can conveniently take into 
account higher order terms. However, we shall restrict ourselves to the first 
two modes for the beam completely off the floor, or to a single touchdown point 
at midspan. 


Solution by Elastica Theory 


For the inflectional case, we have the following equations to describe the 
elastica [2,3]: 


ae as 
p=L == ? (3a) 
B= = Fly) — Fly) » (3b) 
€=2sin-1(k siny) , (3c) 
7 = 7 (EW) — Ell -z (34) 
a ae > cosy = — n = : (3e) 


where s is the arc length (s = 0 at the hinge), (x’, y’) is a coordinate system 
directed along and normal to the resultant load R (Figure 4), ¢ is the angle 


yy 


R (Q-N)/2 
Figure 4 
Coordinates for the elastica. 


made with the x’ axes, & is the modulus of the elliptic functions F’(y) and E(y), 
and y is a parameter, with w= x/2 or — 2/2 corresponding to a positive or 
negative value of € at s = 0. At s = L we have the boundary condition 


¢=-y or ksiny=-—sin (2) (4a, b) 


where y is the angle which the resultant R makes with the horizontal. In 
addition, the x, y coordinates are given by 


x= x'cosy—y'siny, (5a) 


y =x’ siny + y’ cosy. (5b) 
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It is easily verified that € = — y, corresponding to dy/dx = 0, is the only nega 
tive value of ¢ at s =0 for which the constraint condition y(x) 20 can b 
satisfied. This case again is of secondary interest, and we focus attention on 
the case where ¢ > 0 at s = 0. In these cases py = 2/2 and F(yy) = K, E(w) = 
(K and E are complete elliptical integrals). 

This problem is conveniently treated by means of a nomograph. We first 


- 
define dimensionless loads by : 
pw ee G Gu 4 | 

sy aes en ee ee ae Sst 


: 
i 
: 
i | 
ul 
t 


where G = Q — N. Consider 


Es fad : 
idee ire ae Po Ss i Me 8 p | 


By virtue of Equations (4b) and (3b), 9, @2, and @g are all functions of 
sin y. For constant G, y, versus sin y is a radial straight line, while for a constant 
P, @» versus sin y is an ellipse. For a constant k, m3 versus sin y has the ap- 
pearance shown in Figure 5. 

So, for example, given G and P, we first determine the intersection of the 
corresponding y, and @, curves; then we determine the modulus & from the 
(3 curve passing through the intersection, and the angle y from the abscissa 
corresponding to the intersection. Having k and y, the deflected shape is 
completely determined by Equations (3) and (5). In Figure 6 are shown graphs 
of the midspan deflections y(L), computed from Figure 5 and Equation (3e). 

As in the elementary theory, we may identify ‘modes’ of deflected shapes. 
For a value at midspan wy, in the range (7/2) < yz < (3 2/2) there is only the 
inflection point at the support (first mode), for (3 2/2) < pw, < (52/2) there 
is one inflection point at the support and one in the span (second mode), etc. 
On Figure 5, the demarcation between the first two modes occurs at the 
maximum siny value on the yg curves for k < sin (7/4). For k > sin (2/4) it occurs 
at the minimum value of siny lying between the two maxima at siny = 1. At 
the demarcation values the midpoint of the beam is itself an inflection point. 

Figures 5 and 6 allow one to trace the deflected shape as P or Q is held 
fixed and the other varied. For example, for G = 1-0 no deflection is possible 
for P less than Pj, © 1:25, At P = P,,,, the midspan deflection is y(L)/L = 0-4 
and consequently N = 0. Hence G = 1-0 corresponds to Q = 1-0. Increasing 
P above P,;,, one moves along the radial line G = 1-0 on the nomograph 
toward the origin. Two deflections can exist for each P. One corresponds to 
the upper branch of the y(L) curve in Figure 6 (increasing k’ s on the nomo- 
graph), and the other to the lower branch (initially decreasing k's on the nomo- 
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Figure 5 
Nomograph. 
08 | 
| 
t 
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____ Elementary theory 
First mode, 
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-—- Second mote, 
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Figure 6 


Midspan deflection versus P. 
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graph). The upper branch y(L) increases with increasing P until a maximum 
y(L) is reached. (y(L) further decreases because a loop is formed that tightens 


with increasing P.) The lower branch. y(L), however, decreases with increasing” 
P, the behavior being similar to elementary theory. As P further increases, 
this branch changes from a first mode to a second mode configuration (for 

~ sin 3°). Finally at P = 8-24 it leads to y(L) =0 on Figure 6 and the 
y(L) = = 0 curve on the nomograph. Here a reaction force N is mobilized and 


G becomes 0= N. On the nomograph, one may then move along the y( El 
zero line toward the left edge, thus going through a sequence of equilibriu 


states where N goes from zero to Q and P from 8-24 to 8-18. Thus as P is. 


decreased from 8-24 one may either move along the abscissa or the O = 0. 
curve in Figure 6. On the other hand, like elementary theory, elastica theo 


indicates a disjointness in the P for this sequence of equilibrium states. One” 
cannot go continuously along the lower branch Q= = 1-0 curve toa value P> 8-24. | 


However, elastica theory does give upper branch values for P> 8-24. 

Thus, elastica theory differs from elementary theory. In particular, it” 
indicates a bifurcation in the equilibrium states, while elementary theory, also’ 
shown in Figure 6 (dotted lines), does not. The curves typify a snap-through 
effect. That is, one expects that in the lower branch region a sufficient disturb- 
ance will snap the rod through the lower branch configuration into the upper 
branch. 


Stability of the Beam 


The ideas involved in a discussion of the stability of the beam are con- 
veniently illustrated by the finite-degree-of-freedom model shown in Figure 7. 
Torsional springs with spring rate c are placed at the joints within the span to 
provide resisting moments of the form c (8, — 6,). Each rigid bar is of length 
7 BY be 8 

Consider first only deflections that are symmetrical about midspan. These 
are described by two coordinates. For these we take normal coordinates q, 


and qs given by 
: A=ntQ, (6a) 


Gms ea q1 Ee: Lies, (6b) 


The coordinates g, and g, correspond to mode shapes that are analogous to 
sin (7 x/2 L) and sin (3 x x/2 L) for the continuous beam. The origin in a q4, Ye 
plane corresponds to the straight configuration. In terms of g, and q, the 
potential energy V is given by 


V = 0 G2 + 02 93 + Bi G1 — Bod, (7) 


7 
- 
- 


oa 
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The loads P¥ = 5 (3 — V5) ef4, PF = 5 (3 + /'5) c/4 are the analogs of the 
first and third Euler load for a simply supported beam. The constraint of the 
floor requires 6, > 0 and 6, + 0, = 0, or in terms of g, and qo: 


it~, (8a) 


Bid — Bog, 20- (8b) 


We note that at the origin grad V = (f,, — B,) is a vector perpendicular to the 
line f, 7, — Bo 2 = 0 and pointing into the half plane defined by (8b). Hence, 
if one moves away from the origin in a direction such that B, dg, — Bz 4q. > 9, 
that is, so the load Q is raised off the floor, the potential energy increases 
regardless of the magnitude of P. On the other hand, if one moves away from 
the origin along the constraint line f; g, — Bz q, = 0, then 


2 2 a 
ay = Bite BA) ag? = 25 (34/5) [Se- 45] at. 
Thus, if P < P¥ = 25c/4L, the potential energy along the constraint line 
increases, but if P > P¥ it decreases. However, it is easily verified that P¥ is 
just the load that holds the model in equilibrium with the center link on the 
floor and the side links off the floor. That is, P¥ is the analog of the Euler load 
P, for the hinged, built-in column. 


‘h B,9;-Bo4n=9 


9 


V=0 
(<p) 
Q 
Cle, 6) m : m ial oiea 9, +9279 
p i = 
Yj Z (p>p}) —(03'<p<03) 
aes Be 

Figure 7 Figure 8 


Finite-degree-of-freedom model. Geometry of the potential energy. 
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Geometrically, the situation is as shown in Figure 7. The constraint con- 
dition is satisfied in the closed region which is bounded by the lines 6, 9g, — B.q2=90 
and g, + q. = 0 (to the right of the shaded area in Figure 8). The level lines 
V = constant are conic sections in the q,, g, plane. As P increases, the V = 0 
level line is first an ellipse, then a hyperbola, and finally an ellipse again. For 
P < P%, the V = 0 line through the origin lies entirely outside of the region 
satisfying the constraints; for P > P# it is partly inside this region. 

Further, it is readily seen that the above conclusions remain the same if — 
one considers arbitrary deflections satisfying the constraints, rather than — 
symmetric ones. Hence, by the energy criterion, as well as by the ‘neighboring © 
equilibrium’ criterion, one would conclude that the critical load is raised to” 
P¥ by the presence of Q. Likewise, if one assigns a mass m to each rigid bar and | 
assumes conservative impacts with the floor, the total energy of the model is ~ 
conserved. If P < P¥ and the initial deflections and velocities are sufficiently — 
small, the kinetic energy is a positive definite function of 9, 7. and the potential — 
energy is a positive definite function of g,, g.. By the classical kinetic stability — 
criterion, the critical load is thus also P = P#. 

In the case of the continuous beam, we may proced in an exactly analogous © 
way. We write (x) in a set of normal coordinates q,, qo, ..-! 


co 
=> sin 
Gn 2° E; 

cE 


The potential energy of the system is then 


r= 4 Siler (G2) 7] G2 a+ eLesn(h)) 


n=1 


At the origin in the (q,, q, ...) space (straight configuration) 


7 QOLPe _ (nn 
WV = 2) dq, sin (+). (9) 
The constraint at midspan yields 
Ds Gn sin (%") 2.0. (10) 
n=1 ie 


Again, consider first moving from the origin in the (q, qo, .- .) space in a direc- 
tion that lifts the load Q off the floor. In this case, the inequality sign in (10) 
holds and by (9) dV > 0. Next, consider the case where Q remains touching the 
floor. By expanding in the eigenfunctions for a hinged, built-in column, it is 
easy to show that dV > 0 for P below P, (the Euler load for a hinged, built-in 
column), and dV < 0 if P > P,. Hence, by the classical ‘neighboring equili- 
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brium’ and energy criteria the critical load is P,. Likewise, it is readily seen 
that under a sufficiently small disturbance and with conservative impacts, 
oscillations near the straight configuration will remain bounded and small, 
and by the classical kinetic criterion the critical load is P, as well. 

Thus, the classical criteria all agree and, as mentioned at the outset, give 
a result which is physically unsatisfactory. How then is the critical load to be 
defined? From Figures 5 or 6, we see that for an above-the-floor deflection to 
exist at all, P must exceed P,,,,. One can thus be safe if P < P,j,. However, 
except for very large Q, P.in is close to the Euler load P, for a hinged bar, and 
is a very conservative criterion. 

For sufficiently small Q and large P, the elementary theory is statically 
adequate for the lower branches of the y(L) curves in Figure 6, and we may 
investigate the local stability of these configurations by elementary means. 
Such an investigation by the classical methods shows that these configurations 
are locally unstable for P > P,. Thus, in this range of P and Q, we have a 
continuous elastic system that is the analog of a particle which is at the bottom 
of a potential hill, but which is constrained so that it may only move uphill 
(Fig. 9). Under a sufficiently large disturbance, the particle may acquire enough 
energy to go over the top of the hill. For example, from Figure 6 we see that for 


Q = 0-01 and P = 5 the central deflections of the lower branch, locally unstable 
configuration is y(L) = 0-008 L. Hence, for a disturbance of this order of 
magnitude, there is a danger that the beam will snap through the locally 
unstable equilibrium configuration to a large deflection. Further, from Figure 6 
we see that the magnitude of the disturbance required to send the rod into a 
large deflection configuration is arbitrarily small for arbitrarily small Q. In 
terms of the particle analogy, the height h of the hill decreases with decreasing 
Q. Thus, for small Q, one cannot conclude that the critical load is raised to Py, 
unless the expected disturbance is correspondingly small. 


Figure 9 


Single particle analog. 


Likewise as Q is increased for fixed P, the height / of the hill increases. 
Thus, the load Q does in fact have a stabilizing effect. However, it is somewhat 
misleading to say that the effect is simply to raise the critical load as it is 
classically defined. Indeed, we see that a meaningful definition of the critical 
load P in this case must depend not only on the magnitude of Q but on the 


magnitude of the expected disturbance as well. 
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Zusammenfassung 


Das Problem eines horizontalen, gelenkig gelagerten Stabes mit Vertikalkraft Q 
am Mittelpunkt und mit axialer Durchlast P wird behandelt. Eine Durchbiegung 
infolge Q wird jedoch durch eine Bodenebene verhindert. Die gekriimmten Formen 
des Stabes werden durch die lineare Stabtheorie und durch die Theorie der Elastika 
bestimmt. Die letztere ergibt eine Vergabelung der Gleichgewichtslagen, die in der 
linearen Theorie nicht erscheint. Schliesslich wird die Stabilitat des geraden Stabes 
behandelt. Das kinetische Kriterium, das Energiekriterium und das klassische, auf 
nichttrivialen benachbarten Gleichgewichtslagen sich beziehende Kriterium werden 
alle in Betracht gezogen. Alle zeigen, dass eine beliebig kleine Vertikalkraft Q die 
kritische Durchlast P um mehr als das Achtfache erhéht. Das heisst, P geht von 
der Knicklast fiir den gelenkig gelagerten Stab iiber auf die Knicklast fiir den 
halblangen, einerseits eingebauten Stab. Es wird gezeigt, dass die klassischen 
Kriterien in diesem Falle irrefiihrend sind. Nicht nur die Grésse von Q, sondern 
auch die der erwarteten Storung muss beriicksichtigt werden. 
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On the Steady Flow of a Viscous Fluid Past a Sphere 
at Small Reynolds Numbers Using Oseen’s Approximation — 


By Prasuu L. BHaTNaGaR!) and PuRUSHATTAM D. VERMA2), ! 
Bangalore, India ; 


1. Introduction 


The first solutions of the linearized Navier-Stokes equation, which ignored 
the inertia terms for the problem of axi-symmetric flow round spheres, were 
given by SroKEs (1851) [1]*). OszENn (1910) [2] took the inertia terms into ! 
account to a limited extent by assuming that the sphere caused a small pertur- 
bation in the uniform parallel flow and neglected second order perturbation 
velocities. GOLDSTEIN (1929) [3] solved OsEEN’s equations exactly to find a 
series relationship between drag coefficient and Reynolds number R, which 
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compares favourably with experimental results for R < 2, but the predicted 


‘drag is about 10% higher than the experimental value at R = 5. He dit not 


' 


obtain an expression for the stream function from which it would be possible 
to plot the flow pattern. Tomotrka and Aor (1950) [4] discussed the steady 


flow past a sphere at low Reynolds number according to OSEEN’S approximation. 
They discussed the flow patterns around a sphere by considering the flow at 
R=2aU/vy=1, where U is the velocity of the undisturbed stream, @ is the 
radius of the sphere and » is the kinematic viscosity and discovered the presence 
of a stationary vortex-ring behind the sphere. PEARCEY and Mc Hucu (1955) [7| 
carried out the computations for the velocities at Reynolds numbers 1, 4 and 
10. The conclusion of their study is that there is rudimentary evidence of a 
wake at Reynolds numbers even as low as 1, and a wake is definitely present 
at Reynolds numbers of 4 and 10. They found themselves in disagreement with 
TomotiKa and Aot, as they did not find a stationary vortex-ring behind the 
sphere; from their diagrams there were only the barest of hints that it would 


ultimately appear. They drew the non-dimensional graph of fluid velocity in 


radial and angular directions. Since they did not find the stream function they 
could not present the streamline patterns. Their calculations are based on the 
analytic expressions obtained by GOLDSTEIN. 

A finite difference method was used by Kawacutt (1950) [5] to solve the 
Navier-Stokes equation for the flow round spheres at R = 20. He found that 


reversal of flow does not occur on the surface, so that a standing vortex ring 


does not appear for R = 20. This method is extremely laborious but has been 
replaced by Lister (1953) [6] by a relaxation method to compute the flow 
round spheres at R = 0, 1, 10 and 20. Recently JENSON (1959) [8] has applied 
this method to spheres at k = 5, 10, 20 and 40 using a different coordinate 
system. He has discussed the phenomenon of separation of the forward flow 
and development of a circulating wake. He reported the first formation of a 
wake at R= 17. 

It is of interest to examine at what stage of increasing Reynolds number 
a clearly defined wake gets established and to see if the vortex is present. The 
main object of the present paper lies in determining the correct general express- 
ion for the Stokes type stream function for the purpose of computing the flow 
patterns around the sphere. For convenience the essential parts of GOLDSTEIN’S 
analysis necessary for our investigation are first reproduced as briefly as pos- 
sible and then the expression for the Stokes stream function will be obtained. 
The detailed computations have been carried out for the flow pattern around 
a uniformly translated sphere at Reynolds numbers, 1, 4 and 10. This choice 
of Reynolds numbers is made in order to compare the results obtained pre- 
viously by the different authors and to use the values of the coefficients A 
and B,, already obtained precisely by PEARCEY and McHucu on IBM type 
punch-card machines using GOLDSTEIN'S analytical expressions. 


4 
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We find that the disagreement of PEARCEY and McHucu with TomMoTika 


and Aol is probably not due to the truncation of the series at » = 1 but may 
be due to the incorrect Stokes stream functions obtained by Tomotika and Aor 
by integrating only the radial component of velocity. This expression for the 


stream function does not satisfy the angular component of the velocity. A 


general expression for the Stokes stream function has been obtained and on 
detailed computation it has been found that there is no separation of the flow 
at R=1 and R= 4, while at R = 10 the flow is separated and the vortices 
are attached. It is expected that these vortices will be detached as the Reynolds 
number increases. 


2. The General Expression for the Stokes Stream Function 


Let ox, oy, oz be rectangular coordinate axes having the origin O at the 
centre of the sphere which is at rest in a steadily running stream of 
incompressible viscous fluid of infinite extent. We assume that the fluid at 
infinity flows with constant velocity U in the positive direction of the axis of x. 
Then the components of fluid velocity at any point may be written as U + u, v, 
w where wu, v, w are the components of the velocity of perturbation produced 
by the obstacle which become vanishingly small everywhere at a great distance 
from the body. In the absence of external forces the Navier-Stokes equations 
of motion are 


d 0 
or =uV'q—eU 1 Pp, 


q is the velocity vector with components wu, v, w. Neglecting the inertia term 
we get the linearized equations of motion of OSEEN: 


0) >) 
oUt =uV?q-Vp, (1) 
with this we combine the equation of continuity 
V-q=0. (2) 


Let us write U = 2 ky, and it is found that OsEEN’s hydrodynamical equations 
are satisfied by 


it __ 0d TE aye 
doh Mics 2h Ozaa he 

a4 __ od i Oy 
Oe 89k Wet Bes Oops? (3) 
gaz — od ~f Ss 4 


Oz Dik me Ose 
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and 
& p=ou , (4) 
provided that 
ox 
g ee ao ae (5) 
an 
Vb = 0. (6) 


3 Introducing the spherical polar coordinates (7,6, w), the appropriate 
general solution of equation (5) for the case of sphere of radius a, is given by 


y= Ue SB, ya) Pe). (7) 
n=0 


where the B,’s are constants of integration and 


E=kr, w=cosb, 
and 


walt) = (20 +1) (Fe) Knsrnléd (8) 


where K,,, 1/2 (§) denotes the modified Bessel function and P,(u) denotes the 
Legendre function. 
The general solution of equation (6) is given by 


pou yy ee 0) 
n=0 


yn+ fi. 


where the A/’s are constants of integration. 
Taking v,, vg, and v,, the velocity components in the direction of 7, 9 and w 


increasing, we get 


od 1 Ox ¥ COS 0 
% PY) j ) k Oy) . : 
sin 0 fod 1 OX . \ 10 
6 v Y < \ Om 2 ES Ou x J ( ) 


and v,, = 0. 
In order to compute the flow patterns around the sphere we shall now find 


the general expression for the stream function y which is connected with the 
components of velocity by the relations, 


ee 11 
ae ye om ’ ( ) 
and 
: ee (12) 
6 ysin@ or’ 


where yp is the Stokes stream function. 
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sequently the values of y and ¢ from (7) and (9), we have 


av _ _y Sr Aan +) K" Pal) 


Om n-0 5" 

7 *2 aS Hee SB, {%,A6) ah Xn(&)} P, (4) , (13) 
and 
ow cy A! K” P/(u) 
DS ee Rea 

U sin? @ ene B P'(u P 14) 

Sarat wa és i Pilu) =€ Bw} ge). 

Now, 


r) 0 
dy = aa du + ae dé. (15) 


Substituting the corresponding expressions for dy/du and dw/d& from (13) and © 
(14) in (15) and integrating, we find the general expression for yp: 


y7_ (+1) Ay 
y= > Seca yn P,-1() ee eT) ; 
n=0 
U lee) n x 
2p 0-menoSanrin(S AH) se, a9] 
n= r=0 : 
where ¢ is a constant of integration and ¢’s are given by the following set of ; 
equations P,(y) 
$n,0(U) = aia 
ee ee (m+ 1)! P,(u) 
#n,1(4) — 1 Ml [ Paw) (n ie, 1)! 2 , 
i — _ | ™+1)1 Pie) (n+ 2)! P,(u) 
( 7) bn, 2(M) + $1,1(4) - 1)! 2 2! (m — 2)! 2 |: 


UT H) Pale) +2 bali) = Pica a ae Shan a | | 


( 
OF + hal [a A — fi 


PCS eR Se. 19) 8, 8) 6 wwe) te, 6) 8, @ 9 6) fe 
DD) Os UR 9 Oa Re, MOLL, “Bie, Ce ce 6), & 6 W, te>%e [e) ‘alta a “qcini ite’ me! wie auibie meee Ren eae 


(1 — 1) by, nl) + (1 — 1) by na (H) = Foon ! Sil H) oe Pall) is 


and 
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The constant of integration c, can be evaluated if we denote by y = 0 the 
particular stream line which coincides with the axis of symmetry w= + i 
The corresponding stream function for the steady flow past a fixed sphere 
in an infinitely extended fluid can easily be obtained by adding to the right- 
hand side the stream function for the uniform flow, which is given by 


z y= U r* sin?6 , 
See ft 2 
Fi al OS a 


Thus, the required general expression for the stream function is 


1 <1 (n+ 1) A, 
e=-5 UP 0-w)+ Le Sie P,_1(#) — Praale)} 
” cae (1 — p2) eS 0-#) yi w+ 1) B; | area) +c. (18) 
n=0 7r=0 


If we write y, = y/U a, 7, =17/a, nd&=R y,/4, the general expression for 
the non-dimensional stream function y, becomes 


oe) A 


ma tnd-w) +d Saar (Pal) — Pasale} 


et mh (1 — p) eee Oy (2n +1) B, (s ctr) +¢, (19) 


where the Reynolds number R = U d/y is determined by 


3. The Coefficients A, and B, 
The boundary conditions are 
U, = —Ucos§, and vg= Using, at r=a. 


On applying the boundary conditions we get the Goldstein equations for 
finding the coefficients 4, and B,,: 


(tM An 4 TY Ba {xnlEe) Pam(Ee) — mE) Yoon(E0} 
m=0 
me (nm ="1) 


o (21) 
| 0 =0, 2,3,...) 
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“A Lo (ead) 22) 


Eliminating A’’s from equations (21) and (22), we have 


re) ee (23) 
oe a | 0 @>1) 


where 


An, m(So) ad Am(So) Yn, m(§o) %y XmlEo) Pn on (50) 
g Met Dalle fx (Ee) = SP, led em 


0 


sta ll a ah ly 


The detailed computations are carried out by PEARcEY and McHueu for 
determining the values of A, and B, on IBM-type punch card machines for ; 
R= 1,4 and 10. They retained 10 significant figures, and extending to m = 8, © 
m=7for R=1,ton=9,m=8 for R= 4, and to » = 10, m = 9 in the case® 
R= 10. The values of A, and B, obtained by them are used with necessary — 
modifications in the computation of the flow patterns for the Reynolds num- ~ 
bers 1, 4 and 10. i 


Table 2. Values of y 


a % | 1-02 | 1-04 | 1:06 1:08 | 1:10 | 1-20 

0-5 + 2.94996 + 2-96694 
0-6 + 2-13874 + 2-15837 
0:7 +1-20913 +1-22909 
0:8 +0:81875 + 0:81986 +0-82129 + 0-82296 + 0:82710 +0-83925 
0-82 + 0:°74027 + 0:-74128 + 0-74256 + 0-74412 + 0:-74788 +0-75890 
0:84 +0:66104 | +0-66192 + 0:66305 + 0:66444 +0:-66785 +0:67777 
0:86 +0:58101 +0-58180 + 0-58280 + 0-58404 + 0:58709 + 0:59584 
0-88 + 0-50026 + 0-50093 +0:50181 + 0:50286 + 0:-50552 +0:51310 
0:90 + 0:41907 + 0:41963 + 0:42035 + 0-42123 + 0:42320 + 0:42956 
0:92 + 0:33676 + 0:33823 + 0:33782 + 0-33852 + 9-34014 +0:34526 
0:94. +0:25372 + 0:25406 + 0°25448 + 0-25506 + 0-25626 +0-26012 
0:96 + 0:16992 + 0:17009 +0:17041 +0:17077 +0:°17161 +0:17421 
0-98 + 0-08528 +0-:08549 +0-08558 +0-08583 + 0-08619 + 0:08748 
0-99 + 0:04274 + 0:04280 + 0:04288 + 0:04299 +0-04302 + 0:04378 


0-995 +0-02251 | +0-02256 | +0-02259 | +0-02264 | +0-02296 | +0-02306 
0-999 +0-00542 | +0-00542 | +0-00542 | +0-00542 | +0-00544 | +0-00551 
0-9995 +0-00270 | +0-00270 | +0-00270°>| +0-00270 | +0-00270 | +0-00275 
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4. Computation of the Flow Patterns 


Keeping in view the aim of the present investigation, detailed numerical 
calculations have been carried out using the general expression of Stokes 
stream function (19) in the range w=1 to w=0-5 andr=1 to 1, =2 for 
R=1, 4 and 10. The tables of results (2 to 4) give the values of yp to five 
decimal places. They have been rounded off to two decimal places and plotted 
on a graph. Finally the contours of equal values of y are drawn. 


Table 1 
Values of the coefficients A, and B, for k = 1,4 and 10 


R iL 4 10 
n A, Be 28 Be Bs Bo 
(0) + 3-52627 | +0-29563 41-21748 | +2:77817 | +0:71690 | +64-05506 
i +0-45078 | —0-00504 + 0-37926 | —0-45035 | +0-32761 | —28-48695 
2 +0-00222 | +0-00002 +0-01519 | +0-02555 | +0-03363 | + 9:-49657 
3 —0-00002 | —0-03712%10-§ | —0-00132 | —0-00007 | —0-00402 | — 0-85221 
+ +0-00009 | +0-00001 | +0-00022 | + 0-07154 
5 —0-00001 | — 0-00396 
| 6 + 0:00016 | 
for fo = 1 


1-80 1:90 2-00 


42-99752| +3-04059| + 3-09512 4+. 3-16051| + 3-23624 | +3-32183 | + 3-41699} + 3:52151 
+ 2:18941 | +2:23094 | +2:28214| + 2-34247 +2-41149| + 2-48884 | + 2:57422 | + 2:66752 
+ 1-25805| +1:29524 | +1-34006| +1-39209 +1-45092 | +1-51633 | + 1-58807 | + 1:66603 
+ 0:86134| +0-88560 | + 0-91883 + 0-95345 | + 0-99471 | + 1:04058 | + 1:09265 +1-14549 
+0°77898 | +0:80104 | + 0-83120 + 0-86270| +.0-90180| + 0-94183 | +0-98917 +1:03719 
+ 0-69578| +0-71554 | +0:74263 + 077088 | + 080598 | + 0-84190| + 0-88438 + 0:92746 
+ 0-61173 | +0-62920 | + 0-65311 + 0-67804| + 0:70908 | + 0-74078 | + 0:77831 + 0:81638 
+ 0-52685 | +0:54196 | + 0:56265 +.0:58421| + 0-61103| + 0-63849| + 0-67097 + 0°70387 
+ 0-44116] +0:45385 | + 0:47123 + 0:48937 | + 0-51197 | +0-53500| + 0:56233 + 0:59000 
+0:35459| + 9-36482 | + 0-37888 +. 0:39353| + 0:41174| + 0-43036| + 0:45241 + 0-47473 
4+ 0-26717 | +0-27494] +0-28559 +.0:29667 | + 0:31047 | +0-32455| + 0:34123 +0-35814 
+ 0-17898 | +0-18418 | + 0-19132 + 0-19877 | + 0:20803 | + 0-21757 | + 0-22869 +0-24011 
+ 0-08986 | + 0-09247 | + 0:09620 +.0-09994| + 0:10458 | +0-10934] + 0-11498 + 0-12069 
+ 0-04489 | + 0-04630| + 0-04801 +0-04999 | + 0-05224| + 0:05474 + 0:05747 | + 0:06045 
+ 0:02363 | +0-02433 | +0:02519 +0-02617 | + 0:02732| +0-02855 + 0:02994 | + 0-03142 
+0-00562 | +0:00577 +0-:00594| + 0-00614| + 0-00637 + 0:00662 | + 0:00689 | + 0:00719 
+ 0-00284 | +0-00288 | + 0:00294 +0-00306 | + 0-00321 | + 0-00330 + 0:00342 | + 0-00359 
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Table 3. Values of » 

fi hi 1-02 1-04 | 1:06 | 1-08 | 1:10 1-20 
0:5 +1-07049 | +1-10561 
0-6 0-88085 0-91132 | 
0-7 0-68113 0-70587 
0:8 +0:45559 +0:45789 + 0-46040 + 0-46314 0-46786 0-48401 | 
0:82 0-41167 0:41384 0-41609 0:41863 0-42301 0-43778 : 
0-84 0:36738 0-36929 0-37141 0-37368 0-37760 0-39099 iz 
0-86 0-32258 0-32431 0-32618 0-32832 0-33180 0-34370 i 
0:88 0-27733 0:27887 0-28056 0-28235 0:28555 0:29598 | 
0:90 0-23174 0:23318 0-23456 0-23621 0-23873 0-24767 iF 
0:92 0-18591 0-18697 0-18829 0-18950 0-19176 0-19880 | 
0:94. 0-13971 0-14067 0:14166 0-14245 0-14414 0-14978 | | 
0-96 0-09350 0-09391 0:09475 0-09522 0-09653 0-10011 | 
0:98 0:04674 0:04737 0-04723 0-04787 0-04796 0-05061 | 
0-99 0:02424 0-02455 0-02488 0-02521 0-02554 0-02584 | 

| 0-995 + 0-01053 +0:01063 +0-01069 +0-01079 + 0:01090 +0-01141 | 

Table 4. Values of py 

. 4 | 1-02 | 1:04 | 1-06 | 1:08 | 1-10 | 1:20 
0:5 | — 0-55939 — 0-47958 
0-6 — 0-71008 — 0-63509 
0-7 — 0:81539 — 0-75187 
0:8 — 0:84987 — 0:84006 — 0:83013 — 0:82088 — 0:80965 — 0-:76748 
0-82 — 0:82392 — 0°81524 — 0-80635 — 0:79767 — 0-78762 — 0-75031 
0-84 — 0:78817 — 0:78012 — 0-:77233 — 0-76473 — 0-75591 — 0-72309 
0-86 — 0:74147 — 0:73463 —0-72801 —0:72155 — 0-:71404 — 0:68566 
0:88 — 0:68356 — 0-67795 — 0-67254 — 0:66728 — 0-66112 — 0:63791 
0:90 — 0-61112 — 0:60670 — 0:60246 — 0:59831 — 0°59347 —0-57501 
0-92 — 0:52471 — 0:52133 — 0-51809 — 0-51492 —0-51117 — 0:49831 
0-94 — 0:43499 — 0:-41783 —0-41781 — 0:41508 — 0-41195 — 0:40336 
0:96 — 0:29980 — 0:29839 — 0:29710 —0-29579 — 0:29426 — 0-28919 
0:98 — 0:16506 — 0:16445 — 0:16391 — 0-16340 — 0-16272 — 0:16067 
0-99 — 0:07559 — 0:07535 — 0:07514 — 0:07490 — 0:07473 — 0:07366 
0-995 — 0:05849 — 0:05797 — 0:05746 — 0:05695 — 0:05648 — 0-05406 


5. Discussion 


The streamlines are sketched in Figures 1 to 3 for R = 1, 4 and 10. The 
most striking feature of these patterns of flow is that at R = 1 and 4, there 
is no separation of the flow. The flow is approximately symmetrical. The 
streamlines show the change in flow pattern with increasing Reynolds number. 
The actual changes in appearance are very slight. This is in disagreement with 
the results reported by Tomotrxa and Aor (1950) [4]. PEarcEy and McHucu 
(1955) [7] reported a strongly marked wake even at the lowest Reynolds 
number, namely & = 1. The stream lines in Figure 3 shows a separation of the 
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for R = 4 


@ 1-30 | 1-40 | 1:50 | 1:60 1:80 | 1-90 


+1-15356| +1-21347| +1-28461]| +1-36646 | + 1-45857 | + 1-56063 | + 1-67233 | + 1-:79349 


0-95232 1-00308 1:06316} 1-13211 1:20960| 1-29541 1-38934) 1-49114 
0-73855| 0:°77857| 0-82582| 0-87982] 0-94051 1-00768| 1:08117] 1-16076 
0-50887} 0-53523| 0-56991 0:60565| 0:64777| 0-69424| 0-74688| 0-80015 
0-46049| 0-48448] 0:51596| 0-54843] 0-58817| 0-62862| 0-67638|) 0-72461 
0-41151] 0-43298| 0-46129] 0-49032| 0-52589] 0-56213) 0-60486) 0-64795 
0-36191| 0-38083} 0-40584| 0-43148| 0-46282| 0-49476| 0-53238] 0-57031 
0-31165| 0-32817| 0-34973| 0-37188| 0-39890| 0-42641] 0-45886| 0-49171 
0-26097| 0-27473| 0-29289| 0-31156] 0-33418| 0-35740| 0-38451} 0:41197 
0-20975| 0-22086| 0-23558| 0-25056| 0-26886| 0-28737| 0-30936| 033130 
0:15803} 0-16653 0:17762| 0-18889] 0-20268} 0-21662} 0-23301} 0-24981 
0-10583 0-11131 0-11891} 0:12625| 0-13568} 0-14521} 0-15637| 0-16725 
0:05336| 0-05620) 0-05989| 0-06367 0:06826| 0-:07291) 0-07842]| 0-08394 
0-02794| 0-02867} 0-03123] 0-03245}| 0-03546| 0-03713| 0-04058) 0:04268 


| +0-01209| + 0-01286| + 0-01375| +0-01478] +0-01588 | + 0-01714| + 0:01845 | + 0-01993 


+ 0-24661 | + 0-38134 
— 0:00789 | — 0-10553 
— 0-25838 | — 0-17121 
— 0:44497 | — 0:39065 
— 0:46634 | — 0-41832 
— 0:47842 | — 0-43781 
— 0:48057 | — 0-44685 
— 0-47135 | — 0-44349 
— 0-44587 | — 0:-42497 
— 0-40522 | — 0-38907 
— 0:34357 | — 0-33286 
— 0:25627 | — 0-25084 
— 0:14856 | — 0-14660 
— 0:06917 | — 0:06830 
— 0-03810} — 0-03572 


+0-00154 | + 0-12001 
— 0-21502 | — 0-11480 
— 0:-41831 | — 0-34075 
— 0:54937 | — 0-49979 
—0-55581 | —0-51401 
— 0°55557| —0-52013 
— 0-54512 | — 0-51531 
— 0°52435 | —0-49891 
— 0-48660 | — 0-46763 
— 0-43459 | — 0:-42207 
— 0:36146 | — 0-35186 
— 0-26654 | — 0-26204 
— 0-15235 | —0-15074 
— 0-07065 | — 0:06994 
— 0:04272 s 0:04.042 


| — 0-39603 | — 0-30716| — 0-21162| — 0-10887 
— 0-55920 | — 0-48055 | — 0-39748 | — 0-30865 
| — 0-68947 | — 0-62595| — 0-56040| — 0-49143 
— 0-72320 | — 0-68327 | — 0-63874 | — 0:59594 
| —0-71088 | — 0-67543 | — 0-63632| — 0-59861 
| —0-68940 | — 0-65855 | — 0-62512| —0-59269 
| — 0-65726 | — 0-63190] — 0-60347 | — 0-57644 
— 0-61363 | — 0-59283] — 0-57030| — 0-54872 
— 0-55667 | — 0-53985| — 0-52315 | — 0:50679 
| — 0-48556 | — 0-47334| — 0-45998 | — 0-44819 
| — 0-39486 | — 0-38542 | — 0-37804| — 0-36896 
— 0-28405 | — 0-28004 | — 0-27542 | — 0-27149 
—0-15861 | — 0-15719}] — 0-15552| — 0-15416 
| — 0-07304 | — 0-07235| — 0-07171 | — 0-07129 
| —0-05172 | — 0-04943 | — 0-04720 | — 0-04496 


flow at R=10 and the presence of attached vortices. The most reliable 
experimental value for the separation of the flow is R = 14 as given by GARNER, 
Jenson and Keegy (1959) [9]. It is expected that the vortex will gain strength 
as the Reynolds number increases and finally it will get detached. JENSON (1959) 
[8] obtained the separation of the flow at Kk = 17, by a relaxation process for 
flow past a sphere within a pipe. However, the present investigation predicts 
the existence of the separation of flow even at R= 10. Critical Reynolds 
numbers for the separation of flow and for the vortices which become detached 


are under investigation. 
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Figure 1 


14 
Calculated stream-lines past a sphere at R=1. 
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Figure 2 
Calculated stream-lines past a sphere at R = 4. 
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Figure 3 
Calculated stream-lines past a sphere at Rk = 10. 
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Sommaire 


On a obtenu une expression générale pour la fonction de courant de STOKES 
pour l’écoulement d’un liquide visqueux autour d’une sphere. Les lignes de courant 
ont été tracés pour les nombres de ReyNotps = 1, 4 et 10. Le décollement du 
courant fluide n’a pas lieu pour les N. R. = 1 et 4, tandis que pour le N. R. = 10 
les tourbillons attachés se forment dans le sillage. 
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Shear Waves in a Semi-infinite Visco-Elastic Medium due to a 
Sudden Twist Applied at a Point on the Plane Boundary : 


By SuspHENDU K. Datta, Calcutta, India?) 


In this paper, the shear waves produced by a sudden twist applied at a point 
on the plane boundary of a semi-infinite visco-elastic medium at ¢ = 0 have been ~ 
studied. The expressions for the displacement on the plane boundary for both small — 
and large distances from the centre of disturbance have been obtained. The medium : 
is supposed to be isotropic and of the Kelvin type. 


Introduction 


Recently several authors have attempted to explain certain seismic observations 
on the supposition that the material constituting our earth is visco-elastic and of 
the Kelvin type (i.e. the mechanical model is a spring and dash-pot arranged in 
parallel). In a previous paper the present author has studied the shear waves due 
to a torque which is an exponential function of time and which is applied on the 
circumference of a circle on the plané boundary of a semi-infinite visco-elastic 
medium of the Kelvin type. Here the problem when the torque is suddenly applied 
at a point of the plane boundary at the time t = 0 has been solved. Results have 
been obtained in parabolic cylinder functions and at large distances, in terms of 
Hermite Polynomials based on e~*/?, In obtaining these results the method of 
COLLINS (1960)*) has been followed. 


1) The Calcutta Technical School. 
2) See the References, page 563. 
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Tol. XII, 1961 
Equations and solution 

We choose cylindrical polar coordinates (r, 6, 2) with the z-axis directed into 
he interior of the medium, the plane boundary being z = 0, with its origin at the 


entre of disturbance. Then the relevant nonvanishing stress components are 
to) ) (8 “ 


) re (a ay 
eGo NEE ae} \ Soy y 


| 
cy a (w+ a city On,” 
where p, w’ are constants and up is the circumferential component of the displace- 
nent vector. Thus wg will satisfy the differential equation 
ay Rs) 0769 1 Oug up 07u69 1 0?uU9 
al See: et see is 
Cees xr) (on aie ov aa a) oo Ole 
: (1) 
eee. 
ta Wo’ e” 


where g is the density of the medium. 

Let us first suppose that the torque of moment N is applied on the circum- 
ference of the circle y = a and then making a tend to zero we shall obtain the 
solution when the torque is applied at the origin r = 0. 


The boundary conditions are 
T,= —P d(r—a) ot), z=90 (2) 
2 = 0 , = oO, t eS 0 
where e 
| P= ona 
Further, : 
Ug = Ot Oe a 10) 
Here 6(z) is the Dirac 5-function. Following Cortins (1960) we define the two-sided 
Laplace transform 
‘4 co 
jis) = [et 10 at, 
—0o 
and henceforth shall use the notation 
| f(s) =f), H)=Ts), 
to denote the transform and its inverse. This modified Laplace transform of (1) and 
(2) gives : 
Ss 076 1 0Uup ue | aS 3 
(1 + =a) ( or? a y or Oz? 2 thal, 2, 8) » (3) 
Zo,(7, 2,8) =—-Pd(r—a), z=0 | (4) 
(0 Me 2=O. 
Using the Hankel transform defined by 
fo, 2) 
(5) 


u(x, 2, #) ae ug(r, 2, t) Jx(%7) a7, 


0 
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(3) and (4) give 


d*t9(x, 2, S) 5 s? on! ws cy 6 : 
oe fe + (1 + s/o) J ey : a 
oo * 
T,(x, 2, S) =—Pafr hen Lea) dr, z=0 (7) 
0 q 
== Q 5 Z=O, | 
where p | 
— H,(0 2,5) =u (1+ 2) 5 tle, 2,9). (8) 


Solving (6) and using (5), (7), (8) we get 
Pa Tile a) 


tox, 2, S) = mw (1 + s/a) (2 + s?/c? (1 + s/w) (0 
- s2 1/2 ) 
xer|— ("+ amegey) 4: 
On inverting (9) we shall obtain at the limit 
* Ns c (1 + s/@,)t? 
t 0 = 
OE en ick (PRE [3+ veel = 
eee Le Sit sare 
<p OPN CTF slot t 
Writing 
Lt ue(r, 0, t) = u(y, 0,2) and x(r, 0, ) = u(r, 0, s) 
a—0 
we have 
: Ns c (1 + s/@,)1/? 
Oh) == 0 
1G OE) 4naucr (1+ s/w)??? [1 + SY 
11 
aie: os (11) 
| PL e+ slots 
Replacing vy w,/c by € and s by s @, 
1 PTE ae OF aoe eee s 
on i. id “4m more G Bare € Seen dartie 
(12) 
piers aa 
ce 


Again replacing s by s—1, 


t NR Lets 1 : 
leo) faust eg ee" en eee 


0 


E 5-112 


We now expand e° in the form 


oo 
& n 

ess Lame ) ss sone 
nN! 


n=0 
and obtain 


= Eine 
N wme-°* 


u(r, 0, t/w,) et = L 
( /0) 4m mcré 
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Inverting the right hand expression of (13) term by term (which is valid in the 
‘present case) (VAN DER Par and BREMMER, 1950) 


J ‘ ni 
ur, 0, —) gt N 0} ye EM(2 7)" 2 &/8¢ 
"54 4mwcré a4 mn! (wt)? 


x P-»(T57) # a aye (a) Joa pels Gal H(t), 


where H(t) is the step function defined by 


Hivead. f2 9 
=0, #<0 


and D_,,(z) is the Weber parabolic cylinder function of negative integral order — n. 
Replacing ¢ by wy ¢ and writing T’ = w 7%, we get finally 


N a3 oH En (2 T)" ge ee ee 
n=0 | 


u(r, 0, t) 


4Anpoeég ZI reat A eet 
(14) 


x Paz) a te Ba ae a TD») Haye | 


This expression (14) is rapidly convergent for small values of § and vanishes 
when IT —> 0. Moreover it increases for some values of T > 0-and then decreases 
to zero when T > oo. 

| We now attempt to obtain the asymptotic series for u(r, 0, t) when both ¢ andr 
are large. To do this wo use the fact that the inversion of the image when s - 0 

' gives the expression for u(r, 0, t) for ¢ tending to infinity. Replacing s by s « in (12), 
we obtain 


2 v2 
——) = N ws & [« + § @)-8/2 4 oe Poon Oe as)-| eo Sesl(tas)" 


u(r, 0, 4m mck 
(15) 


Wy & 


where we shall take « to be equal to /2/é. To get the asymptotic expansion we shall 
assume s to be small so that |«s| <1. Then, 


- * 


: 1 . , ? (2 — 2)! 
Pees) = So + Dan (sg)P—4 ye ay = (—1)"** Seaman — i) (16) 


Therefore, 
E oo 
me s = = i= Ss = lee Se 3 1 
fae aés+ s? Beane 4 a&s+s 2as?S (17) 
co 
writing S = 3"a,(a s)"~* . Using (17) in (15) 
n=3 
t _ Naja ae ge il aon: toe 
u(r, 0, a=] ~ 4npeé | (itas)?® —& 1l+as > \t> co 


ZAMP X1I/36 
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Again replacing ¢ by ¢+ &a, we obtain 


a, 0, t+é& *) aa re E s p(s) + = 4) f + Sar rts : ( 


Wy % 
where 
p(s) = (Las $= 1— ast > gst .n., 
g(s) = (14+ as\-!=1—as+a%s?— 
and : 


ee =14 Yarra =1—-— 2a,a5s° + 2 a? (aj s* — a, s*) — --- 


Using now the relations, 


oe vol e Flt gn eS = i= ee H en( =) en Flt 

Je Vx Qn |/2+1 x y2 

we obtain 
N wh « eos / 

——— n —u?/4 | 

u(r, 0, t) ~ anpoe |S Fa Po ae pa ont e F (18) 


where u = (T — &)/|/E/2, He,(z) is a Hermite polynomial of order n based on e~*” 


and 
co 


Du Pylt)e*" = as pis) ( + Sa ra ) ; 


1 


co 


Diya" On(t) e-** = g(s) (: +e rt) oe 
1 


0 


It can be shown easily that 


Q,(u) = TIES |Ha( 5) + ay Hel =) ARSE S 


The right hand side of (18) is appreciable for small values of u i.e. for large values 
of € and T (T being very nearly equal to &). 

In conclusion, I express my thanks to Prof. B. Sen, D. Sc., F. N.I., for his 
kind help in the preparation of this paper. 
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Zusammenfassung 


Es werden die Transversalwellen untersucht, welche in einem viskoelastischen 
Halbraum unter einem plétzlich an der Grenzebene angebrachten Moment ent- 
stehen. Insbesondere werden die Verschiebungen an der Grenzflache fiir kleine 
und grosse Distanzen vom Stérungszentrum angegeben. Das Medium wird als 
isotroper Kelvin-Stoff angenommen. 


(Received: March 15, 1961.) 


Mechanical Dislocations and Steady State Processes 


By Ross M. Evan-IwanowskI, Syracuse, N.Y., USA) 


N. I. MusKHELISHVILI [1]*) and later H. PorRITSKY [2] showed that if the temper- 
ature distribution T(¥, y) in an elastic cylinder (plane strain problem) satisfies the 
Laplace’s equation, AT(x, vy) = 0, then the stresses for the auxiliary problem 
(T = 0) are the same as for the original problem (T + 0), (cf. [1], Pp. 157-161). Thus, 
it was proved that the steady state of heat flow is a sufficient condition for this 
to be true. 

Below it is demonstrated that the steady state condition is also a necessary 
condition for the auxiliary problem to have the same stresses as the original one. 

The thermal effect (or an equivalent effect, such as Wigner growth effect 

caused by the fast neutron penetrating some materials) results in the following 


stresses (plane strain problem) : 


0 ov 
re ey TH AO H+ 2 n= 3 pps AO 2p 
(1) 


dv Ou 
el ays oe Soy 
Dividing this problem into two: auxiliary problem (T = 0) and thermal problem 
(T + 0) we write 
n=u ut, v=v' +v* (2) 


where primed variables refer to the auxiliary problem, and the starred variables 


1) Syracuse University, Department of Mechanical Engineering. 
2) Numbers in brackets refer to References, page 564. 
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refer to the thermal problem. Substituting (2) into (1) results in the following: 


Ou , Ou’ Out 
Ov’ Ou* 


Ov’ Ou’ du* ou* 
tay = (Ge + ae) +B a tae : 
To obtain an auxiliary problem independent of the thermal problem, the following 


must be satisfied: 


ou* ov* ou* 
hada ig pk 2 AO*+ 2u 8 Or ae 


AG* 4+ 24 


Ox oy Ox Oy 
which leads to 
ou* Ue Bee Ov* (5) 
Ox Oy Oy CL aie 


These are the Cauchy-Riemann conditions for the analyticity of the function 
do(z) =u* +70". 
From either (4) and (5) it also follows that 
ou* ov* v 
Ox oy 2(A+H) 
Hence T(%, y) is the real part of an analytical function ®(z), such that ; 


T(x, 9). (6) 


$e) =ut Fivt = oa folds, 4 
Re[@(z)] = T(x, y) . (8) 

Moreover . 
AT (x, y) = 0. (°) 


This completes the proof. 

Because of (7), the displacements of the auxiliary problem for the multiply 
connected regions can be multivalued functions (cf. [1], p. 164). This leads to the 
notion of (mechanical) dislocations [4]. To state it differently: introducing in a 
multiply connected body proper types of mechanical dislocations, one can simulate 
(mechanically) the effects of body forces satisfying (1) and (9). Moreover, the last 
proof indicates that only and only such effects are possible to simulate by (mechani- 
cal) dislocations. The significance of this fact, although a negative one, lies in the 
following. The simulation of the thermal (or similar) effects in the plane theory 
of elasticity by application of dislocations can be utilized for the photoelastic 
methods of stress determination [3]. The above theorem states that only steady 
state effects satisfying (1), can be simulated by the type of dislocations considered 
here, i. e., small rotations and small relative displacements. 
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lastiques multiplement connexes, Paris, 


i 


Zusammenfassung 


Es wird gezeigt, dass in einem elastischen, zylindrischen K6rper der stationare 
Warmefluss (oder ein vergleichbarer physikalischer Prozess) eine hinreichende und 
notwendige Bedingung dafiir ist, dass sich die Spannungen aus einem Hilfsproblem 
mit fester Temperatur berechnen lassen. Es folgt daraus, dass sich in der Elastizitat 
ebener Probleme nur stationdre Prozesse bei der Anwendung photoelastischer 
Methoden durch mechanische Dislokationen darstellen lassen. 


(Received: April 25, 1961.) 


d 2 
Remarks on the Differential Equation = = yore 


By J. StEKMANN, China Lake, California, USA?) 


According to a remark by J. v. NEUMANN [1],2) the wellknown Blasius equation 


of order 3 
ft" (n) + 2 f(n) f"(n) = © (1) 


must be reducible to one of first order, because it permits the group of transfor- 
mations 


n>n+m» fn) >R fk 1) - (2) 
Putting . 
= oie At ate _ ae - : 
fmet, hae x(t) , ans oh eae (3) 
-v. Neumann obtains the equation 
Aer). (4) 
dx w(2%—y) ’ 
with the initial condition 
y> oo aS ¥>O. (5) 


This differential equation may be classified as an Abelian differential equation of 
the second kind. By suitable substitution it may be subsumed also in the category 
of the Abel differential equations of the first kind. A formal treatment as an exact 
or total differential equation is possible by means of an Eulerian multiplier. 
In order to study the behavior of the integral curves of the differential equation 
(4), an analog investigation was made by the REAC computer (ayn ore le S. Naval 
Ordnance Test Station, China Lake, California. We note first of all that this first- 


1) U.S. Naval Ordnance Test Station, Research Dpt. 
2) Numbers in brackets refer to References, page 567. 
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order differential equation has singularities at the origin, ie., r= 0, y= 0; a 
infinity, i.c., 7—>0co, y—> OO; and at the intersections of the straight lin 
y +*-+2=0 with the straight line 2% — y = 0 and the y-axis ¥ = 0, ie., the 
points + = — 2/3, y = — 4/3 and ¥ = 0, y = — 2. As these points are approached 
both numerator and denominator tend towards zero and when the magnitude o 
the output voltage of the amplifiers representing these functions approaches thi 
noise level of the amplifiers, the solution becomes meaningless. 


1000 


500 

0 1 

: 

: 

: 

: 

; 
-500 
~1000 


Figure 1 
dy pe ye 2) 


Direction field — 
dx x (2x —y) 


The equation was first coded and scaled to yield a solution within the limits 
— 1000 < * < 1000 and — 1000 < y < 1000. Then the programm was rescaled to 
investigate a very small region near the origin: — 2 <¥ < 2, —3 <y <1. In 
order to check the accuracy of the solutions, the equation for the isoclines y’ = const. 
ba fom ee The results are given in Figs. 1 and 2. They represent the 

irection fie escribed by the differential equation (4). The isocli 

dotted and dot-dashed lines. . o NEE Der 

I would like to thank very much Mr. A. Moline of the NOTS Analog Computing 
Group for providing the programs. 
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Figure 2 
dy y(xtyt2) 


Direction field 
dx x (2x —y) 
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Zusammenfassung 


Verlauf der Integralkurven einer in der Grenz- 


Es wird eine Ubersicht iiber den 
tialgleichung gegeben. Die Lésungen wurden 


schichttheorie auftretenden Differen 
mittels eines Analogrechners berechnet. 


(Received: June 9, 1961.) 
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Ein quadratisch konvergentes Verfahren zur Eigenwertbestimmung i 
bei unsymmetrischen Matrizen I | 


1 


Von Heinz RurTisHausER, ETH, Zirich?) 


7 

Von J.C. F. Francis?) ist kiirzlich ein aus der L R-Transformation [1, 2, 3] ab- 
geleitetes Verfahren zur Eigenwertbestimmung angekiindigt worden, das ebenfalls 
auf der Zerlegung der Matrix in 2 Faktoren beruht, die nachher in umgekehrter 
Reihenfolge wieder zusammenmultipliziert werden. Dieses Verfahren hat den — 
Vorteil, dass die genannte Faktorenzerlegung — im Gegensatz zur L R-Transforma- 
tion — auch bei unsymmetrischen Matrizen nicht infolge Division durch Null 
scheitern kann. 


1. Allerdings ist das genannte Verfahren von FRANcis (im folgenden sei es als 
Qk-Transformation bezeichnet) nur linear konvergent; man kann aber die Kon- 
vergenz durch geeignete Koordinatenverschiebungen erheblich beschleunigen, wie 
in der Originalarbeit von FRANCIS gezeigt wurde. Indessen stellt die Bestimmung 
der Koordinatenverschiebungen eine gewisse Willkiir dar und ist in einem voll- 
automatisch ablaufenden Rechenprozess etwas problematisch. 

Man kann aber — und dies ist der Gegenstand der vorliegenden Mitteilung — die 
QR-Transformation so modifizieren, dass ein im Sinne von GRAEFFE*) quadratisch 
konvergenter Rechenprozess resultiert, der in einem Rechenautomaten vollig 
zwangslaufig und sto6rungsfrei ablauft und damit die «narrensichere» Bestimmung 
der Eigenwerte einer beliebigen quadratischen Matrix gestattet. Dies stellt insofern 
einen wesentlichen Fortschritt dar, als man dies von kaum einem der bisher — 
bekannten Verfahren behaupten kann. 


2. Wenn die gegebene Matrix mit A, bezeichnet wird, so kann die QR-Trans- 
formation im Prinzip durch folgende ALGOL-abnliche *) Rechenvorschrift A definiert 
werden: 


for k:=1 step 1 wntil inf do 


begin 
ihe Zerlege A, in 2 Faktoren Q, x R,, von denen Q; eine orthogonale®), R, 
dagegen eine obere Dreiecksmatrix mit positiv reellen Diagonalelementen ®) 
sein soll; 
a Multipliziere R, und Q, zusammen zu A, ,1 = Ry XQ, 


end 
1) Institut fiir angewandte Mathematik. 

*) Invitational Symposium on Matrix Computations, Gatlinburg, Tennessee, April 1961; er- 
scheint im «Computer Journal». 

3) Der Zusatz «im Sinne von GRAEFFE» deutet an, dass das behauptete Konvergenzverhalten 
fiir alle Eigenwerte gleichzeitig gilt. Im Gegensatz dazu ist ein Konvergenzverhalten «im Sinne von 
Newron», das nur ftir einen einzelnen Eigenwert zutrifft. 

*) Es handelt sich um ein ALGOL-Programm, in dem gewisse zusammengesetzte Anweisungen 
zu Pseudoanweisungen zusammengefasst sind, die die auszufiihrenden Operationen in Worten 
beschrieben. Ausserdem sind die Deklarationen weggelassen. 

») Wir betrachten hier nur reelle Matrizen. Diese kénnen durch orthogonale Transformationen 
zwar nicht allgemein auf Dreiecksform [4], aber immer auf Wintner-Murnaghansche Normalform 
transformiert werden [5]. 

*) Ohne wesentliche Einschrankung der Allgemeinheit kann man annehmen, dass die gegebene 
Matrix nichtsingular sei. Andernfalls kann man diesen Zustand durch kleine Stérungen herstellen. 


f 
, 
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Wie FRANcIs beweist, konvergiert’) die so erzeugte unendliche Folge ahnlicher 
“Matrizen unter gewissen Bedingungen gegen eine Wintner-Murnaghansche Normal- 
form, das heisst gegen eine Matrix mit den Eigenschaften 

aay =seOe, fir + > 7+ 1, 

b) 4; 4-14%i41,,=0 fiir alle = 2,3,...,n—1. (1) 

In beiden Fallen kann man fiir hinreichend grosses k die Eigenwerte der Matrix 
A, aus A,, angenahert ablesen: Reelle Eigenwerte erscheinen als Diagonalelemente 
von A,,; die komplexen Eigenwertpaare erhalt man leicht als Eigenwerte gewisser 
zweireihiger Untermatrizen, die langs der Diagonalen von A, auftreten. 


3. Man kann nun in Analogie zu [1], § 12, die QR-Transformation so modifi- 
zieren, dass nur die Matrizen A,, As, A;, Ay, Ay,... der in Abschnitt 2 genannten 
Matrizenfolge A, in der Rechnung mitgefiihrt werden miissen. Man erhalt dann die 
folgende Rechenvorschrift B: 


begin 
1:  Zerlege A, in 2 Faktoren P,xS,, wobei P, eine orthogonale, S, eine obere 
Dreiecksmatrix mit positiven Diagonalelementen sei; 
2: fork: =1, 2xk while true do 
begin 
<i Bilde das Produkt S,x P, und zerlege dasselbe wieder in 2 Faktoren 
PB, x Se wobei P,, eine orthogonale, S,, eine obere Dreiecksmatrix mit 
positiven Diagonalelementen sei; 


7 ee Bilde das Produkt P,,: = P;, x Ps 

5: Bilde das Produkt Sy,: = S$; Sp: 
end 

end 


Fiir die so erzeugten Matrizen P,, P,, Py, Pg... und Sy, Se, S4, Sq) +. gilt alsdann: 
P,=Q1020s---Qe, Se= Re Res Re-a--- Ke R,, sowie Ayi, = Pe* Ai P,, 
wobei die Matrizen A,,,, Oy, Ry (k = 1, 2, 4, 8, 16, ...) dieselben sind, wie sie durch 
die Rechenvorschrift A erzeugt werden. 

Daher lassen sich die Konvergenz-Eigenschaften der Rechenvorschrift B sofort 
aus den Ausfiihrungen von FRANcIs ablesen, wobei man aber zu beriicksichtigen 
hat, dass fiir die Rechenvorschrift B nur die Gréssen zu den Indexwerten k = 25 
berechnet werden miissen. Dies ist insbesondere im Hinblick auf die Formel 


Agi, Put A, Py (2) 


interessant, welche somit aussagt, dass die Rechenvorschrift B schon nach verhalt- 
nismassig kurzer Zeit eine orthogonale Matrix P, liefert, die die gegebene Matrix A, 
mit geniigender Genauigkeit auf Wintner-Murnaghansche Normalform transfor- 
miert und somit die im Titel genannte Aufgabe im Sinne der numerischen Mathe- 
matik lést (siehe jedoch Abschnitt 5). 

4. Besondere Aufmerksamkeit muss der Wiederzerlegung des Produkts S;, P;, 
in 2 Faktoren P, S, geschenkt werden. Wenn wir nun P, in ein Produkt von 
m = n(n — 1)/2 Jacobische Elementardrehungen zerlegt denken: 


jee Cu Nk ae (3) 


7) Es handelt sich hier nicht um eigentliche Konvergenz, sondern nur um Konvergenz der 
Form nach, das heisst, die auch im limk > oo nicht verschwindenden Elemente der Matrizen Ay 
entsprechen den nichtverschwindenden Elementen einer Wintner-Murnaghanschen Normalform. 


wobei (mit einem gewissen, von 7 abhangigen 1): 


il 
iL 
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} 
5 
: 
7 
= 


COSM, SING; « ~ =) s = Zeile 7 


—sing; cOS%; - ..- - - Zeile 4 + 1 


il 


ist, so kann auch die Operation S, P, = P, Ss ; in Einzelschritte aufgelést werden: 


A 


i 
Sy Uses USPS 
SP u, = 0 SP, 


= A 
SOD a SON 


m 


Die U; sind ebenfalls Jacobische Elementardrehungen, die sich von den entspre- 
chenden U; nur durch den Drehwinkel ¢; unterscheiden, welche durch die Forde- 
rung, dass alle si) obere Dreiecksmatrizen mit positiven Diagonalelementen sein 
sollen, im wesentlichen eindeutig bestimmt sind. Man findet 


eta = 
*%C ef Me . (6) 
wenn ¥%, y, 2, beziehungsweise das (i, 7)-, (2,7 + 1)- und (¢ + 1,2 + 1)-Element der 
Matrix si ~1) bezeichnen. 

Schliesslich wird 


ctg >; = 


S,= Sf", P,=U,U,... U, , (7) 
und damit 
: Pyy = Py Py = Py Uy Um (8) 
sowie wegen (2) und (8): 
beh pie vw ra A 
Agyyy =U,” U_os + U; * Anya U, U,... U, - (9) 


Die Gleichungen (8) und (9) besagen aber, dass man die Matrix A, durch eine 
Folge von Jacobischen Elementardrehungen auf Wintner-Murnaghansche Normal- 
form transformieren kann, die Drehwinkel ¢; werden durch die Formel (6) bestimmt. 


5. Leider kann die Konvergenz des genannten Verfahrens durch 2 Umstiande 
erheblich gestért werden: 

a) Es kann ein zu «disorder of latent roots» (siehe [1], § 4.3) analoger Fall ein- 
treten, der sich darin dussert, dass die Diagonalelemente von S;, sich auch im 
lim k > oo nicht der Grésse nach ordnen. Diese Erscheinung kann man leicht ver- 
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hindern, indem man in der Formel (6) einen Maximalwert fiir ctg p; festlegt,, derart, 
‘dass man bei Uberschreiten dieses Maximalwertes mit M anstelle des wahren 
“Werts von ctg m; in diese Formel eingeht. Ein zweckmassiger Wert fiir M diirfte 
etwa 10! sein. 


b) Falls mehr als 2 Eigenwerte der Matrix A, denselben Absolutbetrag besitzen 


(beispielsweise wenn A, orthogonal ist!), so wird A, fiir k > oo nicht auf die Wintner- 
Murnaghansche Normalform, sondern lediglich auf die allgemeinere Form 


aa 


By 


_transformiert, wobei jede der Matrizen B Eigenwerte gleichen Absolutbetrages hat, 
welche auch Eigenwerte von 4, sind. 


Bevor wir fiir diesen Fall eine Nachbildung des Verfahrens von BRODETZKI- 


SMEALt vorschlagen (das mit einem erhéhten Rechenaufwand verbunden ist), sollen 
erst numerische Experimente abgewartet werden, die zeigen sollen, ob nicht schon 

die durch die Rundungsfehler verursachten Stérungen der Eigenwerte ausreichen, 
um das Eintreten dieses Falles zu verhindern. 


[1] 
[2] 


[3] 
[4] 


[5] 
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Summary 


J.C. F. Francis proposes for the determination of eigenvalues of an arbitrary 


matrix the QR-transformation which is a unitary analogue of the L R-transforma- 
tion. In this process the current approximation matrix is always factorized into a 
unitary and a triangular matrix and the 2 factors are multiplied together in reverse 
order; this yields the next approximant. The present paper shows how this process 
may be extended to a process which transforms a given matrix with quadratic 
convergence orthogonally to WINTNER-MURNAGHAN’S normal form. The process is 


formulated in the real domain. 


(Eingegangen: 4. Juli 1961.) 


Potenzieren von verallgemeinerten Potenzreihen 


Pi 
572 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 
mit beliebigem Exponent | 


Von Micnar ZyczKowsk1, Krakéw, Polen *) 

Bei praktischen Anwendungen von Potenzreihen, zum Beispiel beim Lésen ~ 
von Differentialgleichungen, ist es oft ndtig, die Potenzreihen zu potenzieren. Die ~ 
zweite oder die dritte Potenz lasst sich ziemlich leicht finden mit dem Verfahren — 
der unbestimmten Vorzahlen, welches die Koeffizienten der Anfangsglieder der — 
Reihe nacheinander zu berechnen erlaubt. Einige Anfangsglieder der vierten und 
fiinften Potenz kann man in der Arbeit [7]*) finden. Im Fall von anderen Werten ~ 
des Potenzexponenten werden die Rechnungen oft kompliziert. 

Manche Verfasser (Laska [4], ADAMs [1], RysHik und GRADSTEIN [5]) geben — 
die Rekursionsformel an, welche den Koeffizient A; in Entwicklung. 


(5 “) = Saas ad (1) 
j=0 j=0 


bezeichnet. Angenommen, sei ganzzahlig positiv und a, + 0, haben wir 


i 
Ay=at, Ap=s Dd, miti-j) a, A;_,. (2) 


Es lasst sich leicht nachpriifen, dass diese Formeln bei beliebigem reellem Wert n 
giiltig bleiben, was auch das Radizieren von Potenzreihen erméglicht. 

In vorliegender Mitteilung gehen wir weiter und geben unmittelbare Formeln 
fiir die Koeffizienten A, des Resultates von Potenzieren der Potenzreihen von allge- 
meinerem Typus. Wir befassen uns namlich mit Potenzreihen, deren Exponenten 
auch Bruchwerte annehmen k6nnen; solche Reihen (die zum Beispiel die Bessel- 
schen Funktionen von Bruchordnung bestimmen) werden auch manchmal «verall- 
gemeinerte Potenzreihen» genannt (SMIRNOw [6]). Um von der unbequemen Voraus- 
setzung a) + 0 loszukommen, schreiben wir unsere Reihe in der Gestalt 


Hay pie: ages. (3) 


Diesmal beschrankt die Voraussetzung a) + 0 die Allgemeinheit nicht, da man 
immer ym (Zeiger des ersten nicht Null gleichen Gliedes) entsprechend wahlen kann. 
Die Zahl y kann einen beliebigen reellen Wert annehmen; setzen wir y = 2, so fassen 
wir méglichst kompakt die Reihen fiir gerade (« gerade) und ungerade (w ungerade) 
Funktionen, und zwar fallen die Glieder mit nullgleichen Koeffizienten aus. Im 
Falle, wo w oder » Bruchwerte sind, setzen wir grundsatzlich x = 0, obwohl diese 
Voraussetzung nicht immer notig ist. 
Befassen wir uns mit der Funktion 


ind 5 
Talay Pr ms ; (4) 
7-0 


1) Politechnika Krakowska, Lehrstuhl fiir technische Mechanik. 
*) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 576. 
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wo é eine beliebige reelle Zahl ist. Zuerst die Umformung 


fore) \é 0° g ino) 5 
( Dayar) = mal va v) uf wl ay 7) (5) 
= 7=0 7=0 


wat (6) 


wo 


bezeichnet wird. So wird also das Potenzieren der verallgemeinerten Potenzreihe 
zum Potenzieren einer gew6hnlichen Potenzreihe reduziert. 
Bezeichnen wir 


co oo & 
st) = Sav, GQ = ya, ) On (7) 
j=0 j= 


und suchen die Entwicklung dieser letzten Funktion in der Gestalt 


d 


lo, <) 
G(t) =)7 4,4. (8) 
7=0 
Setzen wir ¢ = 0 in (7) und (8), so bekommen wir direkt 
Ao = ay; (9) 
-‘differenziert man die Funktion (8) -mal nach ¢ und setzt ¢ = 0, erhalt man 


_ a0) . a 


ee n! 


Weil G(t) eine zusammengesetzte Funktion (7) ist, reduziert sich das Problem auf 
die Berechnung der n-ten Ableitung der zusammengesetzten F unktion. Die ent- 
sprechende allgemeine Formel wurde von pE Bruno [2] gegeben und zum Beispiel 
von Goursat [3] und Rysuik und GrapsTEIN [5] zitiert. Diese Formel hat die 
Gestalt 


Ss ' d™G r hy ON ip a go”) Rns 
G@™t) = n! eat Se ae ( | s ($7) aS oy ( -) (11) 
(*) Pa, hy! Bag! Bp, dg™ \1! 2! nl 


ns°* 


wo die Summierung alle ganzzahligen nichtnegativen Wurzeln k,, der Gleichung 


heb Dig tO kgp bo + nha (12) 
. betrifft und S die Zahl der Lésungen dieser Gleichung bezeichnet, dagegen 
m = hyg tt Rag t Bgg +--+ + Fins = MS) - (13) 
In unserem Fall G(#) = [g(t)]}*, also 
d™G eh 
Gee = EE Ema yee, (14) 
und nach (7) haben wir 
g(0) =i! a,. (15) 
Daher 
(n) . n! é (E ad, 1) (E —m+ 1) an” ais gts “ ans (16) 
ee Slo! 4 OD hs at ne 
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und endlich 


is q 

=1)...(é—m-+41 a ne ae k 

Ay ci » ae . = ne mite gee ame et 2, 3... on 
sie 2380 "es ns" 


Diese Formel mit der Formel (9) zusammen bilden die gesuchten allgemeinen — : 
Formeln fiir die Koeffizienten der mit beliebigem Exponent potenzierten Potenz- — 
reihe, denn nach (4), (5), (6), (7) und (8) erhalten wir 


co \é ied g | 
i (Sa) Spe sg: (18) 
j=0 7-0 


Weil & eine beliebige reelle Zahl ist, umfasst diese Formel auch das Radizieren — | 
einer Potenzreihe und das Dividieren einer Zahl durch eine Potenzreihe. | 

Das Suchen der ganzzahligen nichtnegativen Wurzeln der Gleichung (12), : 
welches zur Anwendung der Formel (17) nétig ist, ist fiir kleine » sehr einfach. — 
Pir y=) hat man nur A, 1) also S = 1: fur 2-2) hat man fh, — yes ] 
hy = hy = 0, Oder hy = hy. = 0, hy = gg = 1, also S= 2, usw. Die vier ersten — 
Glieder A, haben also die Gestalt : 


A,=4a., 
A,=€a*-1a,, 
2 TY ae + na Se 
Ae Bed a (€ ~ 2) ag-9 gS 4 tise a bg a4 a ab? a, (19) 
Ae a (eo 3) a 2) (€ — 3) af-4 qt + als ats — 2) ai aa, 

“| us aks 0 ay + rs rags © oF aoa Gan | 


Bei héheren Werten fiir x, die ubrigens in Anwendungen seltener auftreten, wachst 
schnell die Zahl der Summanden S in Ausdriicken fiir A,,: beim = 5 hat man S = 7; 
bel v= 6, S=11; bein=7, S= 15. Fiir diese Werte m sind alle Lésungen der 


Gleichung (12), mit entsprechenden Werten m zusammen, in den Tafeln 1, 2 und 3 
zusammengestellt. 


Tafel 7 
Die Wurzeln der Gleichung (12) fiir » = 5 
s ky ky ks ky ks m 
il 0 0 0 0 1 1 
2 it 0 0 1 0 2 
3 0 1 il 0 0 2 
4 2 0 1 0 0 3 
5 1 2 0 0 0 3 
6 3 1 0 0 0 4 
WS 5 0 0 0 0 5 


j 
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| Tafel 2 
F Die Wurzeln der Gleichung (12) fiir = 6 
; Ss ky ky ky ky ks ke m 
1 0) 0 0 0 0 al i 
2 ib 0 0 0 1 0 Ds 
3 10) 1 0 1 0 0 2 
4 0 0 va @) 0 0 2, 
5) 2 0 0 1 0 0 3 
6 iL iL il 0 0 0 3 
? 0 3 0) 0 0 0 3 
8 3 0 1 0 0 0 4 
9 2 2 0 0 0 0 4 
10 4 it 0 0 0 0 5 
esa) 6 0 0 0 0 0 6 
Tafel 3 
Die Wurzeln der Gleichung (12) fiir » = 7 
s hy ky Rs Ry Rs Re ka mH 
1 0 0 0) 0 0 0 1 1 
v4 1 0 0 0 0 ul 0 2 
3 0 il 0 0 u 0 0 2 
4 0 0 iT it 0 0 0 2 
5 2 0 0) 0 1 0 0 3) 
6 iL 1 0 il 0 0 0 3 
7 a 0 2 0 0 0 0 3 
8 0 Z 1 0 0 0 0 3 
9 5 0 0 “ 0 0 0 4 
10 2 a ik 0 0 0 0 4 
Tl uy B | 0 0 0 0 0 4 
12 4 0 al 0 0 0 0 5 
TS Ss 2 0 0 0 0 0 5 
14 ss] 1 0 0 0) 0 0 6 
15 =" 7 0 0 0) 0 0 0 7 


Als Beispiel berechnen wir einige ersten Glieder der Entwicklung der Funktion 


F(x) = [Jurs(*)}?” 
in die verallgemeinerte Potenzreihe. Die Besselsche Funktion Jj),(*) besitzt die 
Entwicklung 


(giana ee ee -) 
Jurs(*) -S rye ar rene rasy \ “16 + 396 35840 


daher « = 1/3, »y = 2, und da wir € = 3/2 angenommen haben, erhalt also die ge- 
suchte Reihe die Gestalt 


co 
[hira(*)}*” = A; yll2+25 
j=0 


Die Formeln (9) und (19) benutzend berechnen wir A, und endlich 


1 9 405 3159 
1/2 5 12 49/2 18 /2 ae -) ; 
y2 (4/3) (« 32°: * 14336 2293760 


Sirs)? = 
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Die Konvergenz dieser Reihe ist gut; zum Beispiel fiir + = 1 nach Behalten de 
vier Glieder der Reihe, weil I” (4/3) = 0,8929795. [9], bekommen wir [J,;3(1)]?/* ~ 
0,624818, anstatt des exakten Wertes (WarTSON [8]) [Ji;3(1)]*/2 = 0,624 836. 
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Le ET AD ECE 


Summary 


The paper gives general formula (17) for coefficients of any real power & of a — 
generalized power series (4). Four first coefficients are given by (19). An example 
of application to Bessel functions is described. i 


(Eingegangen: 10. Mai 1961.) 
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Eighth National Symposium on Reliability and Quality Control 
in Washington, D. C., USA 


In Washington, D.C. (USA) findet vom 9, bis 11. Januar 1962 das 8th National 
Symposium on Reliability and Quality Control statt. Betriebssicherheit und Quali- 
tatstiiberwachung sind Gesichtspunkte, die in heutigen elektronischen Gerdaten, 
welche Tausende und oft Zehntausende von Schaltelementen enthalten, grésste 
Bedeutung erlangt haben. Vom einwandfreien Funktionieren einer elektronischen 
Anlage kann die Produktion eines ganzen Betriebes, ja im Kriegsfall ganzer Landes- 
teile abhangen. Der Zweck dieser Tagung, welche vom Institute of Radio Engineers 
(IRE) und anderen wissenschaftlichen Vereinigungen durchgefiihrt wird, ist, Spe- 
zialisten, welche auf dem Gebiet der Betriebssicherheit und Qualitatsiiberwachung 
uber Erfahrung verfiigen, zusammenzufiihren. In Vortragen und Diskussionen wird 
ein ausgiebiger Gedankenaustausch gepflegt. Programme und Anmeldeformulare 
sind zu beziehen bei A. H. DRAYNER, Publicity Chairman, The Martin Co., Balti- 
more, Md., USA. A, P. SPEISER 
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-_ IUTAM-Symposium tiber Kreiselprobleme 


Die Internationale Union fiir Theoretische und Angewandte Mechanik wird 
vom 20. bis zum 23. August 1962 in Celerina (Oberengadin) ein Internationales 
Symposium tiber Kreiselprobleme (Theoretische Grundlagen, Anwendungen, Kreisel- 
geratetechnik) abhalten. 

Interessenten wollen sich an Prof. Dr. H. ZIEGLER, Eidgendssische Technische 
Hochschule, Ziirich, wenden. 


et 


4 


( 


Die Erstellung eines internationalen mehrsprachigen Verzeichnisses 
von Begriffen und Fachwértern tiber Rechenanlagen 


Durch die Zusammenarbeit von zwei internationalen Organisationen wird zur 
Zeit ein mehrsprachiges Verzeichnis von Begriffen und Fachwortern tuber die 
Terminologie der automatischen Datenverarbeitung erstellt, welchem internationale 
Bedeutung zukommen wird. Die zwei Organisationen, namlich die Internationale 
Féderation fiir Datenverarbeitung (IFIPS) und das Provisorische Internationale 
Rechenzentrum (PICC), sind iibereingekommen, das Verzeichnis nach seiner Voll- 
endung der Internationalen Organisation fiir Standardisierung (ISO) zur Verfiigung 
zu stellen. A. P. SPEISER 


Vereinigung Europdischer IBM-1620-Benutzer 


Es wird die Griindung einer Vereinigung Europaischer IBM- 1620-Benutzer 
geplant, nach dem Vorbild der 1620 Users Group in den USA. Die Hauptaufgabe 
einer solchen Vereinigung wiirde Meinungs- und Informationsaustausch der Benut- 
zer und zukiinftigen Benutzer des IBM-1620-Elektronenrechners sein. Interessenten 
werden gebeten, sich mit Herrn Ing. Dr. Hans Tompa, European Research Associa- 
tes, 95 rue Gatti de Gamond, Briissel 18, Belgien, in Verbindung zu setzen. 


ee 
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Experimental Nuclear Physics. Vol. III. Von E. SEGRE (John Wiley & 
Sons, Inc., New York 1959). 811 Se lOehis.; » 23.—- 

Les physiciens attendaient avec impatience le 3éme et dernier volume de la 
série Experimental Nuclear Physics éditée par E. SEGRE (John Wiley) et dont les 
volumes I et II ont paru il y a plus de cinq ans. C’est maintenant chose faite et 
800 pages nous présentent les aspects expérimentaux et théoriques les plus récents 
des phénoménes de désintégration radioactive. 2 

Le volume est divisé en cinq parties d’inégales longueurs. La premiere (partie 
VIII de la série), d’une cinquantaine de pages, est due a E. SEGRE et traite des 
caractéristiques générales du processus de désintégration. Bien que la matieére en 
soit classique, l’adjonction de tables de standards radioactifs, le traitement dé- 
taillé des phénoménes de fluctuation et un apercgu des méthodes de mesure des 
périodes les plus courtes en fait un résumé utile et fort bienvenu. 
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La seconde partie (partie IX) rédigée par G.C. Hanna présente, en quelqu 
200 pages, un excellent exposé de la radioactivité « sous tous ses aspects. C'est am 
notre connaissance le seul compte-rendu récent de ce domaine qui soit si complet 
et présenté d’une facon aussi systématique. Nous ne doutons pas qu’il sera appré- 
cié des spécialistes comme il le mérite. I] discute d’abord les relations énergétiques 
puis le probléme essentiel des facteurs gouvernant la vie moyenne en relation avec 
les différents modéles nucléaires, enfin il passe en revue les méthodes de détection 
et de spectroscopie « pour se terminer par une table fort compléte des émetteurs 
« et de leurs caractéristiques. 

M.Devutscu et O.KoroED-HaNsEN assument la lourde responsabilité de 
traiter dans les 3éme et 4é¢me parties les phénoménes d’émission y (partie X) et de— 
désintégration B (partie XI). C’est grace a leur connaissance parfaite du sujet et 
a leurs contributions personnelles dans ce domaine qu'il leur est possible de pré-— 
senter de facon quasiment exhaustive et up to date des domaines en aussi rapide — 
évolution. C’est ainsi que la non conservation de la parité et ses conséquences ont _ 
pu étre traitées de fagon satisfaisante en dépit de leur intervention récente. Cepen-_ 
dant l’effet MOssBaverR dont on connait l’importance capitale pour la spectros-_ 
copie nucléaire a haute résolution est trop jeune en date pour avoir pu étre pris en : 
considération. Il semble bien que ce soit, de nos jours, le sort inévitable des traités © 
qui se veulent complets d’étre déja dépassés dans l’une ou l’autre partie au mo-— 
ment méme ow ils paraissent. 

Quoi qu’il en soit, le présent compte-rendu est d’une richesse et d’une clarté 
qui devraient satisfaire les plus exigeants. 

La derniére partie (partie XII de la série) est consacrée aux accélérateurs de 
particules et a pour auteur E.M.McMILLan. On comprend que les éditeurs aient 
désiré inclure dans leur ouvrage cet aspect dominant pour ne pas dire écrasant de 
la physique expérimentale moderne. L’énorme accumulation de types nouveaux 
et de techniques en développement a contraint l’auteur 4 un exposé de nature 
plutét historique et descriptif en dépit de ses 150 pages, qui tranche avec le 
caractére précis et quantitatif des autres articles. Le lecteur en appréciera 
certainement l’ampleur et l’abondante illustration mais ne pourra se retenir de le 
considérer comme une adjonction un peu en marge de l’ouvrage. 

Comme ses prédécesseurs, le volume III de Experimental Nuclear Physics — 
bénéficie d’une présentation de haute qualité et d’une typographie parfaite. Nous 
sommes persuadés qu’il clot dignement une série dont la valeur est indiscutable 
et qui rendra longtemps encore d’inestimables services 4 tous les physiciens, 
théoriciens comme expérimentateurs, engagés dans l’exploration du noyau. 


J. ROSSEL 


Aufsatze und Vortrage tiber Physik und Erkenntnistheorie. Von Wo trF- 
GANG PauLt. (Verlag Vieweg Sohn, Braunschweig 1961). 183 S., 7 Abb.; DM 19.80. 

Dieses tiberaus gewichtige Buch von kaum 200 Seiten enthalt 20 Aufsitze und 
Vortrage, die, bis auf eine Ausnahme, zwar alle schon friiher, zum Teil in englischer 
Sprache, gedruckt vorlagen, aber erst in dieser Form einem breiteren Leserkreis 
zuganglich werden. Das Hauptziel der Veréffentlichung, die vom Verfasser noch 
selbst geplant war, lasst sich am besten durch ein Zitat aus dem Vortrag iiber Die 
philosophische Bedeutung der Idee der Komplementaritdét (S. 34 des Buches) be- 
schreiben. Dort steht als einleitender Satz: «Die Veréffentlichung dieses Vortrages 
erfolgt in der Hoffnung, mit diesem kleinen Beitrag diejenigen grésseren Bestre- 
bungen zu férdern, welche im allgemeinen das Ziel verfolgen, die verschiedenen 
Teildisziplinen, in welche unsere Geistigkeit auseinandergefallen ist, einander 


ey 


. 
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wieder naher zu bringen.» Und wenig spater: «(Aber) heute scheinen mir die Voraus- 

-setzungen fiir ein erneutes Einverstandnis der Physiker und Philosophen tiber die 
‘erkenntnistheoretischen Grundlagen der wissenschaftlichen Naturbeschreibung 
erfiillt zu sein.» 

Und nun folgt man mit anhaltender Spannung dem grossen Selbstgesprach*) 
des Verfassers, der wie kein anderer die tiefgreifenden Wandlungen der Physik in 

“unserer Zeit bewusst und kritisch erlebt hat; iiber die ordnende Kraft der Idee der 

-Komplementaritaét im physikalischen und in ausserphysikalischen Bereichen. 

Als die grossen Gestalter der heutigen Wissenschaft treten NrELs BoHR, WERNER 
HEISENBERG, vor allem aber ALBERT EINSTEIN und in bemerkenswerter Weise 

-C.G. June auf. ARNOLD SOMMERFELD sind zwei Beitrage gewidmet und ein 
menschlich ergreifender Nachruf geht an PAUL EHRENFEST. 

Besonders hervorzuheben ist an allen urspriinglich deutsch geschriebenen Bei- 
tragen der hervorragend schéne sprachliche Ausdruck, der die Ausgewogenheit mit 
der vollkommenen Klarheit verbindet. : 

Eine besondere Stellung nimmt der letzte Aufsatz iiber die altere und neuere 

Geschichte des Neutrinos ein. Er ist ganz der Physik gewidmet und beschreibt vor 
allem aktuelle Experimente iiber den f-Zerfall und daran anschliessende theore- 
tische Uberlegungen. Es war WoLFGANG Pav ti zwei Jahre vor seinem Tode ver- 
génnt, die endgiiltige Bestatigung des von ihm 25 Jahre zuvor postulierten Neu- 
trinos und dessen Rolle bei der Paritatsverletzung noch zu erleben. Dieses Erlebnis 
ruft beim Referenten die Bemerkung Schopenhauers (Aphorismen zur Lebens- 
weisheit, Kapitel IV iiber den Ruhm) in Erinnerung, welche lautet: «Im Alter giebt 
es keinen schéneren Trost, als dass man die ganze Kraft seiner Jugend Werken 
einverleibt hat, die nicht mitaltern.» 

Wir sind dankbar, in diesen Aufsatzen eine einzigartige Quelle zur Entwicklung 
der wichtigsten naturwissenschaftlichen Ideen der ersten Halfte des 20. Jahr- 
hunderts zu besitzen. Res Jost 


Proceedings of the 10th Annuai International Conference on High 
Energy Physics. Rochester, 25. August bis 1. September 1960. Herausgegeben von 
E. C. G. SupaRSHAN, J. H. Trntor und A. C. MELIssiINos (Interscience Publishers, 
New York 1960). 890 S., 621 Fig.; Sao. 50. 

Publier un volume de 890 pages, excellement illustré et d’une typographie 
parfaite deux mois seulement aprés la conférence, constitue un tour de force qu’il 
convient de souligner. 

| Cet ouvrage constitue en quelque sorte une somme des connaissances actuelles 
sur les interactions des particules fondamentales aux hautes énergies. L’organi- 
sation de la matiére refléte celle qui a prévalu lors de la conférence: quatre sessions S 
portant sur les aspects expérimentaux et théoriques de l’interaction forte des pions 
et nucléons et des particules étranges, ainsi que sur les interactions faibles, com- 
posent la premiére partie (pages 3 a 619) avec 130 contributions présentées a la 
conférence et une quarantaine d’articles destinés uniquement aux comptes-rendus. 
Il n’est pas possible dans une bréve revue de citer les auteurs de ces communications 
parmi lesquels se trouvent les meilleurs spécialistes du moment. 

La seconde partie (pages 623 4 753) qui englobe les quatre sessions R reprend 
les mémes problémes mais sous forme d’exposés d’ensemble présentés par des rap- 
porteurs particuliérement qualifiés pour dégager clairement Vétat actuel des 
expériences les plus relevantes de la théorie des interactions fortes et faibles, plus 


1) Die dem Verfasser eigentiimliche, eindriickliche Form des Vortrages war dem Selbstgesprach 
viel naher verwandt als der Auseinandersetzung mit einem Partner (dem Horerkreis). 
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dentes. 

Les problémes de la vie moyenne du pion neutres, de la photoproduction et de 
la diffusion x — p au voisinage des résonances supérieures ainsi que la conservation 
du spin isobarique, critére de l’indépendance de charges des forces nucléaires, ~ 
sont résumés par J. ASKIN; en particulier l’expérience semble confirmer toujours 
mieux le dernier point. Une autre question importante est l’interaction 2 — 7 sur 
laquelle les expériences actuelles ne peuvent naturellement fournir que des indica-~ 
tions indirectes et ont beaucoup de peine a apporter des précisions valables. | 
G. C. Wick présente une courte revue de ]’état actuel de l’application des théories | 
de dispersion en insistant sur la complication intrinséque de ces méthodes. Les © 
données sur les phénoménes intervenant entre nucléons et antinucléons sont dis- | 
cutées par M. Lévy dans leurs aspects généraux, et par K. SyMANZIK en ce qui © 
concerne les procédés mathématiques utilisés pour l’analyse des processus d’inter- : 
action. La revue des données experimentales sur l’interaction forte des particules _ 
étranges est assurée par M. ScHwarrtz et D. MILLER et celle des aspects théoriques © : 
par P. T. Matruews. L’état actuel du probléme des interactions faibles est pré- _ 
senté sous les trois aspects de leur intervention entre les particules normales — 
(V.L. TELEGp1) entre les particules étranges (M. SCHWARTZ) et sous l’angle — 
théorique (A. L. GOLDBERGER et L. B. Okun). La difficulté de séparer proprement 
les interactions faibles et fortes et leurs relations trés complexes sont particulié- — 
rement soulignées ainsi que le caractére encore trés problématique du boson inter- — 
médiaire invoqué pour expliquer l’interaction faible. 

_ La troisiéme et derni¢re partie du volume est consacrée au compte-rendu des 
3 sessions P de la conférence. La session P, traite de la structure des particules — 
«élémentaires» étudiée par diffusion cohérente des électrons de haute énergie 
(facteur de forme du proton et du neutron) par collision entre faisceaux d’électrons | 
par des mesures de moments magnétiques et enfin par diffusion des muons de 
grande énergie (expérience du CEern). Il s’agit d’un ensemble de 13 communications. 

La session P, (10 articles) présente de nouveaux résultats dans le domaine des 
énergies ultra-élevées, en particulier les travaux du CERN avec le synchrotron a 
proton de 25 GeV et des recherches par les rayons cosmiques (Moscou, Bristol, 
Chicago). 

Une dizaine de communications consacrées 4 la théorie des particules élémen- 
taires constitue la session P,; on y reléve avec intérét l’utilisation dans ce domaine 
également des concepts théoriques de la superconductivité montrant que leur 
efficacité ne s’étend pas seulement au domaine du solide et de la structure du 
noyau. 

Tel qu’il nous est présenté le volume des Proceedings contitue un ouvrage de 
référence que tout physicien, méme sans étre spécialiste des hautes énergies, peut 
et doit consulter avec un trés grand profit. J. RossEL 


. 


Statistics Manual. Von L. Crow, A. Davis und W. MaxrFieLp, Research 
Department, U.S. Naval Ordnance Test Station (Dover Publ. Inc., New York 
1960). $1.55. 

In diesem Handbuch sind saémtliche heute bekannten statistischen Tests, die 
in der Forschung benutzt werden, sehr klar und iibersichtlich zusammengestellt. 
Fiir jeden Text werden die entsprechenden Begriffe erklart und seine Anwendung 
anhand von Beispielen gezeigt. Zahlreiche Tabellen und Diagramme erleichtern die 
Benutzung der Tests. Das Buch kann allen Wissenschaftern, die Statistiken, For- 
schungsresultate usw. analysieren miissen, warm empfohlen werden. W. SAXER 
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der Mechanik und Physik 


Herausgegeben von Prof. Dr. Pu. FraNxK, Boston, und .Prof. Dr. 
R. von MisEs, Cambridge/USA, unter Mitarbeit namhafter Wissen- 
schaftler. 


Nachdruck der 2., vermehrten Auflage. 1961. DINA5. 


Leinen DM 66.-. Beide Bande werden nur geschlossen geliefert. ~ 


Die Differential- und Integralgleichungen 


Erster, mathematischer Teil: XXIII, 916 Seiten mit 84 Abbildungen. 
Zweiter, physikalischer Teil: XXIV, 1106 Seiten mit 110 Abbildungen. ; 


Dieses von der Fachwelt seit langem erwartete Werk erscheint nunmehr als Nachdruck der zweiten, 
vermehrten Auflage zu einem dusserst niedrigen Preis. Der erste Band bringt in erster Linie eine Dar- 
stellung aller wesentlichen in der Physik bendtigten mathematischen Methoden, wahrend der zweite 
Band eine Beschreibung ihrer Anwendung auf die verschiedensten Probleme der Physik gibt. 


Messen und Rechnen 
in der Physik 


Von Prof. Dr. U. St1tLz, Braunschweig. 2., verbesserte und erwei- 
terte Auflage. 19,5 x 27,5 cm. XII, 472 Seiten mit 6 Abbildungen, 
55 Tabellen und 35 Tafeln. 1961. Leinen DM 78—. . 


Die Neuauflage ist gegeniiber der ersten wesentlich iiberarbeitet und erweitert. Dieses wird besonders - 
deutlich im Kapitel «Ergdnzungen». Das Literaturverzeichnis wurde von 950 auf 1350 Nummern erwei- - 


tert, wabrend das Sachregister jetzt 2000 gegeniiber 1700 Anfiihrungen enthalt. 
Inhalt: Einfiihrung und Begriffsbestimmung / Mechanik / Warme und Strahlung / Elektrizitat und 


Magnetismus / Akustik und Phonometrie — optische Strahlung und Photometrie / Werte fiir Konstanten /— 


Erginzungen: Grossen und Grdssenarten; Gréssen gleicher Art und Grdssen gleicher Dimension; Ver- 
haltnisgréssen und Verhdltniseinheiten — Internationales Einheitensystem; «Metrisches System»; ver- 


einheitlichtes yard und pound — Ebener und riumlicher Winkel als Verhaltnisgréssen — Vom Prototyp- 


meter zum Wellenlangenmeter — Entwicklung der Zeitintervalleinheit Sekunde; Weltzeit, Ephemeriden- 
zeit, sAtomzeit» — Trag 


e und schwere Masse — Technische Krafteinheit — Die Druckeinheiten Bar und _ 


Torr — Neudefinition der internationalen Tafel-Kalorie und der British thermal unit; Bezugszustand der — 


internationalen Dampftafeln — Temperatur und Temperaturskalen; Gold- und Platinpunkt; photo- 
metrisches Strahlungsaquivalent — Wandlungen des Molbegriffes; «Teilchenmenge» als Grundgréssenart 
und als Verhaltnisgrésse - Zusammenflihrung der chemischen und der physikalischen Atomgewichts- 
skala zu einer vereinheitlichten relativen Atommassenskala / Tafeln. 
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